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1 Introduction 

The theory of heavy-flavour production covers in reality a wide variety of phe- 
nomena, which are at times only loosely related from a theoretical point of view. 
In ref. [Nason92] an ample discussion was given of the next-to-leading order (NLO) 
correction to the hadroproduction of heavy flavours for the total and single-inclusive 
cross section and for the double-differential distributions. Both hadroproduction and 
photoproduction have been considered in the literature, and NLO calculations of 
heavy-flavour production in deep-inelastic scattering have been available for a long 
time. Some very special topics, such as the multiplicity of heavy flavours in jets, have 
also been discussed in the previous volume. 

In this chapter we will discuss and review recent progress that took place since 
the publication of ref. [Nason92] . Many theoretical and experimental issues on heavy- 
flavour production have been clarified, and new problems have arisen. First of all, 
the top quark has been observed, and its cross section turns out to have roughly the 
predicted magnitude. In the future, both at the Tevatron and at the LHC, it will 
be possible to study more details of the ti production mechanism. From the point 
of view of perturbative QCD this is a very interesting possibility, since theoretical 
predictions for the top are very reliable. This is in fact the only case in heavy-flavour 
production when radiative corrections to the total cross section are indeed small. 
They are of the order of 15% at the Tevatron energy. The first studies on kinematical 
distributions of the top have already appeared, showing a qualitative agreement with 
QCD predictions. We will review in this work the present status of the theoretical 
calculation of top cross section and distributions. The importance of these studies is 
twofold. On the one hand one would like to test the underlying strong interaction 
dynamics, and on the other hand, deviations from QCD predictions may hint at the 
effects of physics beyond the Standard Model. 
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The study of bottom production at the Tevatron has undergone a considerable 
improvement, thanks to the introduction of vertexing techniques. The inclusive 
transverse-momentum distribution of h mesons has been a long-standing problem, 
since the measured cross section used to be several times higher than the theoreti- 
cal calculation. Experimental refinements of the measurement have reduced the cross 
section by a considerable amount. Nevertheless, the data still tend to lie on the upper 
side of the theoretical prediction. Although it is clear by now that the hard QCD pro- 
duction mechanism is basically the correct one, it appears that some details remain 
to be understood. Vertexing techniques have also allowed additional measurements 
of correlations between the h and the h quark, adding evidence for the validity of the 
hard production mechanism. 

Fixed-target studies of heavy-flavour production have provided a wealth of data 
on charm. Total cross sections, single-inclusive distributions, correlations between the 
quark and the antiquark have been measured in both hadro- and photoproduction. 
The theoretical apparatus of perturbative QCD is in this case at its very limit of 
applicability, because the charm mass is very close to typical hadronic scales. Thus, 
effects of non-perturbative origin will very likely play an important role. Conversely, 
it is hoped that these effects may be better understood by studying charm production. 
Although modern fixed-target experiments have considerably improved the situation, 
many open problems remain in this field. All experimental results are in qualitative 
agreement with perturbative QCD calculations, thus supporting the "hard" nature 
of charm-production phenomena. However, several quantitative deviations from pure 
QCD are observed. It is interesting to see whether simple models of non-perturbative 
phenomena, such as fragmentation effects, intrinsic transverse momenta, and dragging 
effects, may be sufficient to explain these deviations. We will discuss these problems 
at length. As of now, in our opinion, several different scenarios are possible, and more 
experimental data may help to discriminate between them. 

The electron-proton collider HERA has begun to produce data on charm photo- 
production and electroproduction. In this case, thanks to the wide kinematical range 
potentially available, the possibility exists of performing QCD studies on the proton 
and photon structure functions. 

Bottom and charm production studies at LEP have reached a remarkable stage of 
precision. Accurate measurements are now available for the h production rate and for 
its fragmentation function. Furthermore, measurements of the gluon splitting rate 
into h and c pairs have been performed. 

On the theoretical side, new developments have taken place, in relation to the 
experimental progress that has been achieved. Thus, top cross section calculations 
have been revisited, and the problem of threshold effects has been re-examined and 
considerably clarified. The resummation of perturbative effects at large transverse 
momentum has been investigated, because of its relevance to the interpretation of 
the CDF and DO data. With a similar motivation, production of heavy-fiavoured jets 
has been considered. Lastly, the computation of heavy-flavour electroproduction has 
been extended to next-to-leading order, in view of comparisons with HERA data. 

We have limited the scope of the present work to open heavy-flavour production. 
The problem of quarkonium production is of a rather different nature, not directly 
related to the topics touched in our review. We have tried to be as complete as 
possible. However, we have skipped topics that are, in our opinion, less important 
for the purpose of testing the underlying production dynamics. 
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2 Fixed-target production 

Heavy-flavour production has been studied extensively in fixed-target experiments, 
with both hadron and photon beams. The typical centre-of-mass energy is in the range 
10-40 GeV, where the bottom cross section is rather small. Therefore most of the 
available data are on charmed-hadron production. 

In this section we will present a comparison between existing fixed-target data 
for the hadroproduction and the photoproduction of heavy flavours and NLO QCD 
predictions. Because of the small value of the charm-quark mass, the perturbative 
expansion may not be reliable, due to higher-order terms of the perturbative expansion 
and to non-perturbative effects. We will discuss these problems and the importance 
of such contributions. More details are given in refs. [Mangano93] and [Frixione94a] . 

2.1 Total cross sections 

A list of experimental results on total cross sections for D-meson hadroproduction 
is presented in table |l[ We used these results to estimate the total charm-pair cross 
section acs- In most cases fixed-target experiments give their cross sections with the 
Feynman-x cut Xp > 0. The ratio crcc/o'cc{x'p > 0) has been evaluated theoretically 
[Mangano93] and is approximately equal to 1.6 in pion-nucleon collisions and to 2 in 
proton-nucleon collisions; it is nearly independent of the heavy-quark mass and beam 
energy (at least for rric between 1.2 and 1.8 GeV, and beam energy Eb between 100 and 
1000 GeV). Therefore, the total DD cross section is obtained from a (D/D, Xp > j 
by dividing by 2 to get the pair cross section from the single-inclusive one, and 
multiplying by 1.6 (for pion) or 2 (for proton) to account for the partial coverage of 
the x^-range. 

The contribution to ace from Ac and Ds production must also be included. We 
use (see ref. [Aoki92] and references therein) 



a{D, 



a{D^ + D+) 

a(Ae) 
a{D^ + D+) 



0.2, (2.1) 

0.3 (2.2) 



(here particle means also antiparticle). Measurements of the cross sections for Dg 
and Ac have recently been presented by the E769 [Alves96] and WA92 [Adamovich96] 
collaborations, and values consistent with eqs. (|2.1|) and (|2.2| ) have been found. There- 
fore, to obtain ace from the total cross section for DD production, we have to multiply 
by a factor of 1.5. 

The WA75 collaboration [Aoki92] does not present single-inclusive cross sections 
for D mesons, but quotes directly the cc cross section 

ace = 23.1±1.3+^:i|/ib, (2.3) 
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(GeV) 


a (total) or a^ {xp > 0) 




D+/D^ 


pN 

[KodamaQl] 


800 


a{D^/D^) = 38 ±3 ±13 
a{D+/D-) = 38±9±14 


38 ±10 


1.0 ±0.6 


pN 

[BarlagSS] 


200 


(T+{D/D) = 1.5 d= 0.7 ±0.1 


1.5 ±0.7 





PP 

[AmmarSS] 


800 


^(^7^0) = 22+? ± 5 
a{D+/D-) = 26 ±4 ±6 


24 ±6 


1.2 ±0.6 


PP 

[AguilarSS] 


400 


a{D^/D^) = 18.3 ±2.5 
a{D+/D-) = 11.9 ±1.5 


15.1 ±1.5 


0.7 ±0.1 


pN 

[Alvcs96] 


250 


^+(^7^0) = 5.6 ±1.3 ±1.5 
a+{D+/D-) = 3.2 ± 0.4 ± 0.3 


8.8 ±1.5 


0.57 ±0.22 


Tl-p 

[AguilarS5a] 


360 


a+{D^/D^) = 10.1 ±2.2 
(r+{D+/D-) = 5.7 ±1.6 


12.6 ±2.2 


0.6 ±0.2 


R-N 

[Barlag91] 


230 


a+iD^D^) = 6.3 ±0.3 ±1.2 
a+{D+/D-) = 3.2 ± 0.2 ± 0.7 


7.6 ±1.1 


0.5 ±0.2 


R-N 

[BarlagSS] 


200 


^+(^7^0) =3.4l°i±0.3 
a+(D+//^-) = 1.7l°J±0.1 


A 1 +0-6 


0.5 ±0.1 


R-N 

[Kodama92] 


600 


a+{D^/D^) = 22.05 ± 1.37 ± 4.82 
a+{D+/D-) = 8.66 ± 0.46 ± 1.96 


24.6 ±4.3 


0.4 ±0.1 


R-N 

[Alvcs96] 


210 


a+{D^/D^) = 6.3 ±0.9 ±0.3 
a+{D+/D-) = 1.7 ±0.3 ±0.1 


6.4 ±0.8 


0.27 ±0.06 


R-N 
[Alvcs96] 


250 


a+{DyD^) = 8.2 ±0.7 ±0.5 
a+{D+/D-) = 3.6 ± 0.2 ± 0.2 


9.4 ±0.7 


0.44 ±0.06 


R + N 
[Alvcs96] 


250 


a+{DyD^) = 5.7 ±0.8 ±0.4 
a+{D+/D-) = 2.6 ± 0.3 ± 0.2 


6.6 ±0.8 


0.46 ±0.09 


R-N 
[Adamovich96] 


350 


a+(D75°) = 7.78 ± 0.14 ± 0.52 
a+{D+/D-) = 3.28 ± 0.08 ± 0.29 


8.8 ±0.5 


0.42 ±0.05 



Table 1: Experimental results on total cross sections for charm production. The 
notation (7+ represents the inclusive cross section for positive rapidity. The pair 
cross section a^i) includes our corrections for possible Xp cuts, as discussed 
in the text. The experimental results have not been corrected for the updated 
D^K-R branching ratios [Barnett96]. 
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for a 350 GeV vr" beam on emulsion, assuming an atomic mass dependence A°'*^ 
(this result was not inserted in table 0). The E789 collaboration [Leitch94], using an 
800 GeV proton beam colliding on a Be or Au target, studied the nuclear dependence 
of the D-meson cross section, finding an A°' behaviour, with a = 1.02 ± 0.03 ± 0.02. 
As a by-product, they obtain 



a 



[d^/D^) = 17.7 ± 0.9 ± 3.4 /ib. (2.4) 



An analogous measurement for the charged-D cross section is not expected^. 

From table ^ we see that most of the experimental collaborations report separately 
the total cross sections for charged and neutral D-meson production. In the last 
column of table ^ we show the ratio R between these two values. For pion-nucleon 
collisions, the agreement between various collaborations is fairly good (except for the 
result of ref. [Alves96] obtained with a beam energy of 210 GeV), and the errors 
on R are moderate; the D^ cross section is found to be about twice as large as the 
D^ cross section. A simple model for estimating the charged-to-neutral D cross 
section ratio is the following. One assumes isospin invariance in the c^D and c^D* 
transition. Furthermore, one assumes that the D cross section is one third of the D* 
cross section, due to the counting of polarization states. Using then the published 
values of the D*-^D branching ratios [Barnett96], the result is roughly 

Observe that this number has changed with respect to the value of 0.43 quoted in 
ref. [Frixione94] , because of changes in the measured branching ratios. This number 
is compatible with the values reported in table |1| for vrA^ collisions, especially for 
the most recent data. For pN collisions, the same simple argument should hold. 
Experimental measurements do not give a unique indication; some results seem to 
give comparable total cross sections for charged and neutral D. We do not find any 
reasonable explanation of this fact. For example, the HERWIG Monte Carlo program 
[MarchesiniSS, Marchesini92] gives roughly the same R value for the proton and for 
the pion beams. On the other hand, the experimental data for the proton beam are 
much less clear than in the case of n'N collisions, the uncertainty on R being quite 
large. The two R values with the smallest uncertainty, namely those reported in 
ref. [Aguilar88] and in ref. [Alves96] , indicate a behaviour similar to the one observed 
in TT^N collisions. 

In fig. |l| we plot the cc and bb cross sections, computed in QCD at NLO, as func- 
tions of the beam energy, for tt^N collisions. The same quantities are shown in fig. ^ 
for a proton beam. The cross sections are calculated using the parton distribution sets 
of ref. [Martin95a] for the nucleon, and the central set SMRS2 [Sutton92] for the pion. 
The default values of the charm and bottom masses are 1.5 and 4.75 GeV respectively, 

^M. Schub, private communication. 
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Figure 1: Pair cross sections for b and c production imr N collisions 
versus experimental results. 



and the default choices for the factorization scale fip and the renormalization scale 
Aifl are 

lip = 2mc, Hr = nic (2.6) 

for charm and 

fiF = ^R = rrib (2.7) 

for bottom. 

The bands in the figures are obtained as follows. We varied fiR between half 
the central value and twice this value. The factorization scale fip was also varied 
between mh/2 and 2mh in the case of bottom, while it was kept fixed at 2mc in 
the case of charm. This is because the adopted parametrizations of parton densities 
are given for Q^ larger than 5 GeV^. The bands shown in the figures are therefore 
only an underestimate of the uncertainties involved in the computation of charm 
production cross sections. We verified that considering independent variations for 
the factorization and renormalization scales does not lead to a wider range in the 
bottom cross section for the energies shown in the figures. We also show the effect of 
varying rric between 1.2 GeV and 1.8 GeV, and mi, between 4.5 and 5 GeV. 

The proton parton densities of ref. [Martin95a] are available for a wide range of 



Aqcd values, corresponding to a^ 



iniy 



values between 0.105 and 0.130. The bands 



shown in figs. ^ and |] for bottom production are obtained by letting Aqcd vary in 
this range. In the case of charm, values of Aqcd corresponding to as{mz) above 
0.115 induce values of as{rnc) too large to be used in a perturbative expansion. For 
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Figure 2: Pair cross sections for b and c production in pN collisions 
versus experimental results. 



this reason, the upper bounds on charm production cross sections are obtained with 
as{rnz) = 0.115. We point out that, by varying Aqcd, one is forced to neglect 
the correlation between Aqcd and the pion parton densities, which were fitted in 
ref. [Sutton92] with Af^ = 122 MeV. 

Experimental results on bottom production at fixed target have been reported in 
refs. [BariQl, BasileSl, Kodama93, Catanesi89, Bordalo88, Jesik95] for pion-nucleon 
collisions. Recently, the first measurement of the bottom cross section with a proton 
beam has become available [Jansen95]. Preliminary results for pN collisions have 
also been presented in ref. [Spiegel96]. These results are shown in figs. |l| and |^. We 
made no effort to correct the data in order to get the 66-quark cross section when the 
i?-hadron cross section was reported. 

As can be seen, experimental results on total charm cross sections are in reasonable 
agreement with theoretical expectations, if the large theoretical uncertainties are 
taken into proper account. We can see that the hadroproduction data are compatible 
with a value of 1.5 GeV for the charm-quark mass. In the case of bottom production, 
the spread of the experimental data is almost as large as that of the theoretical 
predictions. Both currently available data points lie in the theoretical band. 

We remind the reader that many puzzling ISR results in pp collisions at 62 GeV 
cannot be currently explained (see the review in ref. [TavernierS?]), in particular the 
large A^ production rates reported in refs. [Bari91, BasileSl]. 

Total cross sections for charm production have also been measured in photopro- 
duction experiments. In fig. H the relevant experimental results of refs. [Alvarez93, 
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Figure 3: Pair cross sections for c production in jN collisions versus 
experimental results. 



Anjos89, Anjos90, Bellini94] are shown in comparison with NLO QCD predictions. 
As can be seen, the theoretical uncertainties are smaller in this case than in the 
hadroproduction case. Again, a charm mass of 1.5 GeV is compatible with photopro- 
duction data. It should be stressed that some of the experimental results are totally 
incompatible with one another. Until these discrepancies are resolved, it will not be 
possible to use the data to constrain physical parameters. For example, while the 
E687 data are inconsistent with a charm mass of 1.8 GeV, this mass value cannot be 
excluded because of the E691 data. 



2.2 Single-inclusive distributions 

Many hadroproduction and photoproduction experiments have measured single- 
inclusive Xp and Pt distributions for charmed hadrons in ttN collisions [Kodama92, 
Barlag88, Barlag91, Aguilar85a, Alves92, Adamovich92, Aoki88, Aoki92a] and in pN 
collisions [Kodama91, Barlag88, Ammar88, Aguilar88, Adamovich92]. Distributions 
are expected to be more affected by non-perturbative phenomena than total cross 
sections. For example, an intrinsic transverse momentum of the incoming partons, 
and the hadronization of the produced charm quarks, may play an important role in 
this case. We will therefore try to assess the impact of such phenomena by means of 
simple models. 

Recent measurements of single-inclusive differential cross sections have been per- 
formed at CERN by the WA92 collaboration [Adamovich96] , which uses a tt" beam 
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Figure 4: The single-inclusive p'^ distribution measured by WA92 (left) 
and E769 (right), compared to the NLO QCD predictions, with and with- 
out the inclusion of non-perturbative effects. 



15 



of 350 GeV colliding with isosinglet nuclei, and at FNAL by the E769 collabora- 
tion [Alves96a] with pion, proton and kaon beams of 250 GeV on isosinglet targets. 

In fig. § we show the comparison between the single-inclusive p"^ distributions 
measured by the WA92 and E769 collaborations in ttN collisions and the theoretical 
predictions^ The solid curves represent the pure NLO QCD predictions for charm 
quarks. The dashed curves show the effect of adding to the perturbative results an 
intrinsic transverse momentum of the incoming partons {kj. kick). This procedure is, 
to a large extent, arbitrary. We implemented it in the following way. We call Pt{QQ) 
the total transverse momentum of the pair. For each event, in the longitudinal centre- 
of-mass frame of the heavy-quark pair, we boost the QQ system to rest. We then 
perform a second transverse boost, which gives the pair a transverse momentum equal 
topr((5Q)+/cT(l)+^T(2); fc-r(l) and/i;r(2) are the transverse momenta of the incoming 
partons, which are chosen randomly, with their moduli distributed according to 



1 dN 
Nd^ 



{kl) 



exp{—kl/{kl)). 



(2.8) 



Alternatively, one may proceed as in the previous case, but giving the additional 
transverse momentum /ct(1) + A;t(2) to the whole final-state system (at the NLO in 
QCD, this means the QQ pair plus a light parton), and not to the QQ pair only. We 
verified that the two methods give very similar results for (A;J) smaller than about 
2 GeV2. 

Another non-perturbative effect that must be accounted for is the hadroniza- 
tion process. This effect can be described by convoluting the partonic cross section 



^In this subsection and in the foUowing one, the relative normahzation of experimental distribu- 
tions and theoretical curves has been fixed in order to give the same total rate. 
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with a fragmentation function, for which we choose the parametrization proposed in 
ref. [Peterson83]: 

D(x) = ^, (2.9) 

^ ' x[l-l/x-e/(l-x)]' ^ ' 

where N is fixed by the condition / D{x)dx = 1. We used e = 0.06, which is the central 
value quoted in ref. [Chrin87] for D mesons. The fragmentation process degrades the 
parent charm-quark momentum, and thus softens the p^ distribution. 

From fig. ^, we see that the effect of the kj. kick on the predictions for bare quarks 
results in a hardening of the p"^ spectrum, overshooting the data. On the other hand, 
by combining the kj, kick with (/c^) = 1 GeV^ and the Peterson fragmentation, the 
theoretical predictions undershoot the data (dot-dashed curves), although no serious 
inconsistency can be inferred given the size of the experimental uncertainty. Taking 
the data at face value, one can notice that larger values for {k'^) and for nic would 
improve the agreement. We have checked that, in order to reproduce the WA92 and 
E769 data, an average intrinsic transverse momentum (fc^) = 2 GeV^ is needed. This 
value is exceedingly large with respect to the typical hadronic scale of a few hundred 
MeV. On the other hand, the data are well described (dotted plots) by the theoretical 
curves including fragmentation and kj. kick with (fc^) = 1 GeV^, if the larger value 
of rric = 1.8 GeV is adopted. 

We conclude that, when using the central value of the charm mass, rric = 1.5 GeV, 
favoured by the total cross section measurements, the theoretical results for the p"^ 
spectrum can describe the data well, only if a large k-^ kick is applied to the fragmented 
curve. If a larger value of the charm mass is adopted, a more moderate and physically 
acceptable fc^ kick is enough to get a good agreement with the measurements. There- 
fore, there seems to be a potential discrepancy between theory and experiments in 
the Py spectrum in charm hadroproduction. From a different point of view, however, 
the discrepancy may be interpreted as the signal that some of the other theoretical 
assumptions are not totally sound. For example, the Peterson fragmentation function 
may not be suitable to describe the hadronization process in hadronic collisions; the 
data would suggest a function more peaked towards the x ~ 1 region. Moreover, 
higher-order perturbative corrections may also play a role, especially in the low-pj. 
region. 

We now turn to the Xp distribution. The experimental measurements of refs. 
[Kodama92, Barlag88, Barlag91, Aguilar85a, Alves92, Adamovich92, Aoki88, Aoki92a] 
(for ttN collisions) and of refs. [Kodama91, Barlag88, Ammar88, Aguilar88, Adamovich92] 
(for pN collisions) find on average a behaviour that is harder than the perturbative 
QCD result for bare quarks (see ref. [Frixione94] for details) . As in the case of the p"^ 
distribution, some description of the hadronization phenomena should be added to 
the perturbative calculation in order to compare it with the data. This problem was 
considered in ref. [Mangano93], where the hadronization phenomena were studied us- 
ing the parton shower Monte Carlo HERWIG. In ref. [Mangano93] the conclusion was 
reached that the combined effects of perturbative higher orders and non-perturbative 
(partonic intrinsic transverse momentum and hadronization) contributions eventually 
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Figure 5: Experimental Xp distributions for D mesons, compared to 
the NLO QCD prediction for cfiarm quarks. 



result in a hardening of the Xp distribution for bare quarks. There it was also argued 
that the usual approach of complementing the perturbative calculation with a frag- 
mentation function in order to describe the Xp distribution is completely unjustified, 
since the factorization theorem only holds in the large-p^ region. 

The most recent experimental results of WA92 [Adamovich96] and E769 [Alves96a] 
for the D-meson Xp distribution in iiN collisions are, instead, in agreement with the 
perturbative QCD distributions for bare quarks, as can be seen from fig. ^ This 
is also consistent with previous findings of E769 [Alves92]. The agreement is quite 
satisfactory in almost the whole range considered for both experiments (which are 
performed at different beam energies). The shape of the theoretical curve shows a 
mild sensitivity with respect to the value of the charm mass (all curves have been nor- 
malized to the data). Recent results for a proton beam are also available [Alves96a], 
and an agreement with theoretical predictions similar to the one displayed in fig. ^ is 
found. 

Almost all of the experimental collaborations have observed, in pion-nucleon col- 
lisions, the so-called leading-particle effect, that is, an enhanced production at large 
Xp of those D mesons whose light valence quark is of the same flavour as one of the 
valence quarks of the incoming pion. In ref. [Frixione94] it was shown that the QCD 
predictions are in better agreement with the data for non-leading particles than with 
the data for the full D-meson sample which, as mentioned before, displays on aver- 
age a harder behaviour than the theory. The difference between the leading and the 
non-leading sample may be an indication that non-perturbative phenomena (such as 

^The pion beam used by the E769 collaboration [Alves96, Alves96a] is a mixture of vr^ (70%) 
and 7r+ (30%). To produce the data appearing in fig. ^, which contain the contributions of both Z?+ 
and D^ , the distributions obtained with the different beams have been combined together, since no 
statistically significant discrepancy has been found between them [Alves96a]. 
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Figure 6: Experimental p^ distribution compared to the NLO 

QCD predictions (rric = 1.5 GeV), with and without the inclusion of 
non-perturhative effects, in 'jN collisions. 



colour-drag effects) are present in the production of leading particles. Recent experi- 
mental results on the asymmetry between D~ and D~^ production were presented in 
ref. [Aitala96], for tt~N collisions with Ei, = 500 GeV. It was found that perturbative 
QCD, the instrinsic charm hypotesis [Combridge79, BrodskySO], and the standard 
Monte Carlo simulations cannot describe the measured quantities. The data can 
only be acceptably reproduced with a special tuning of Monte Carlo models imple- 
menting beam-dragging effects. Further studies of colour dragging and hadronization 
phenomena have been presented in refs. [Aitala96a, Carter96]. 

In the past, the Xp and p^ distributions have been fitted with the functions 



da 



Ail-x^r; 



da 



dp. 



Ce-''P' 



(2.10) 



and the data have been presented in the form of measured values for the parameters 
n and b. A possible way of comparing the data with theoretical predictions would 
be to fit theoretical distributions in the same way, and then compare the values of 
the parameters n and b. This was done in ref. [Frixione94] . However, as pointed 
out there, this procedure is not satisfactory. In the case of the p^ distribution, for 
example, it was found that the function 



da 
dp^ 



C 



bml + Pj 



(2.11) 



provides an excellent fit to the theoretical curve. Since, as discussed before, the data 
show the same qualitative behaviour as the theoretical predictions over the whole p^ 
range explored, the form of eq. (|2.11|) is clearly preferable, as was explicitly shown 
more recently in ref. [Alves96a]. 
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Single-inclusive distributions for charm production have also been measured in 
photon-nucleon collision experiments. In the case of photoproduction, we expect 
QCD predictions to be more reliable than in the hadroproduction case, since only 
one hadron is present in the initial state (see refs. [Ellis89, Smith92, Frixione94a] for 
a detailed discussion). In fig. ^ we show the p^ distributions measured by the E687 
[Bellini94] and by the E691 [Anjos89, Anjos90] collaborations. We also show the NLO 
QCD prediction for bare quarks, and the QCD prediction supplemented with Peterson 
fragmentation and an intrinsic transverse momentum for the incoming partons with 
different values of (fc^). It is interesting to notice that, in this case, the fragmentation 
effect, combined with a moderate intrinsic transverse momentum of the initial-state 
partons, is sufficient to reproduce the experimental distributions. Contrary to what 
happens in the hadroproduction case, the p"^ distribution is now less sensitive to the 
choice of the (fc^)- To show this fact clearly, we have also presented the prediction 
for (kl) = 2 GeV^. 

The E687 collaboration recently performed a measurement [Frabetti96] of the 
asymmetries between charm and anticharm states in photon-nucleon collisions, find- 
ing no compelling evidence for Df and A^ asymmetries, and sizeable asymmetries 
for D^, D^ and D*~^ states. The very tiny asymmetries predicted by NLO QCD 
cannot account for these results, which are most likely to originated from some non- 
perturbative effects, as discussed in the case of hadroproduction. It was found that, 
by properly tuning Monte Carlo programs based on the Lund string fragmentation 
model, the data could be reproduced. 

In summary, a remarkable amount of experimental information on the single- 
inclusive distributions for charm quarks is currently available, both in hadroproduc- 
tion and in photoproduction. Further data with increased statistics are expected in 
the near future [Kaplan94]. Some modelling of non-perturbative effects is needed 
in order to describe the transverse-momentum spectrum of charmed hadrons. The 
inclusion of a fragmentation function tends to give too soft a Pt spectrum. Assuming 
the presence of a primordial transverse momentum of the colliding partons, the p^ 
spectrum becomes harder, and a better (albeit not fully satisfactory) agreement with 
the data is found. Longitudinal (xp) distributions are measured to be harder than 
the theoretical predictions. Similarly, asymmetries in the longitudinal distributions 
between charm and anticharm states are too large to be explained by perturbative 
QCD alone. Simple models of hadronization, including colour-dragging effects, can 
give a satisfactory description of the data. 

We remind the reader that the E653 collaboration presented in ref. [Kodama93] 
the first measurement of single-inclusive distributions for bottom production at fixed 
target (the azimuthal correlation of the bb pair was also studied). The data appear to 
be in qualitative agreement with QCD predictions. More investigations with improved 
statistics will be welcome in this area. 
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2.3 Double-differential distributions 

Many experimental results on correlations between charmed particles in hadro- 
and photoproduction have been obtained by different experiments (see for exam- 
ple refs. [Aguilar85b, Adamovich87, Aguilar88, BarlagQla, KodamaQla, Aoki92a, 
Alvarez92, Frabetti93]); these reported distributions of the azimuthal distance be- 
tween the charmed hadrons, the rapidity difference, the invariant mass and the trans- 
verse momentum of the pair. A detailed comparison of these results with QCD 
predictions is performed in ref. [Frixione94] . More recently, new measurements of the 
azimuthal distance and pair transverse momentum for charmed mesons have been 
presented by the WA92 collaboration [Adamovich96a] . In what follows we will fo- 
cus on the distribution of A0, defined as the angle between the projections of the 
momenta of the pair onto the transverse plane, and of the transverse momentum of 
the pair priQQ)- We will discuss whether NLO QCD predictions can describe the 
available experimental data. 

In leading-order QCD the heavy-quark pair is produced in the back-to-back con- 
figuration, corresponding to A0 = tt and PriQQ) = 0. NLO corrections, as well as 
non-perturbative effects, can cause a broadening of these distributions, as illustrated 
in refs. [Mangano93] and [Frixione94a] . 

We have chosen, as an illustration for hadropro duct ion, the cases of the WA75 
and the WA92 results, which have both been obtained in tt'N collisions at the same 
energy. Eh = 350 GeV. Let us first consider the A0 distribution. In fig. ^ we show 
(solid curves) the NLO result superimposed on the data of two experiments. The 
charm mass was set to its default value, rric = 1.5 GeV. In both cases, we see that 
the experimental data favour a much broader distribution than the pure NLO QCD 
result for charm quarks. 

One should, however, take into account also non-perturbative effects, as in the 
case of single-inclusive distributions. We have computed the A0 distribution in per- 
turbative QCD with an intrinsic transverse momentum of the incoming partons as 
described in subsection ^^ (the use of a fragmentation function has no effect on the 
A(f) distribution, since it does not affect momentum directions). The dashed and dot- 
ted curves in fig. |^ correspond to the NLO prediction, supplemented with the effect of 
an intrinsic transverse momentum with (A;^) = 0.5 GeV^ and (fcj) = 1 GeV^, respec- 
tively. We see that with {k"^) = 0.5 GeV^ it is impossible to describe the WA75 and 
WA92 data. This conclusion differs from the one of ref. [Frixione94] for the WA92 
result. This is because in ref. [Adamovich96a] the WA92 collaboration has improved 
the study of correlations with respect to ref. [Adamovich95] by considering a wider 
set of correlation variables and by improving the statistics by a factor of 5. WA92 
and WA75 data now appear to be consistent with each other. As is apparent from 
fig. 0, the acceptable value of {k"^) = 1 GeV^ is required to describe the data. 

We point out that this result has to be interpreted in conjunction with the adopted 
values of the input parameters entering the calculation. We normally choose the 
following defaults: 
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Figure 7: Azimuthal correlation for charm production in irN collisions: 
NLO calculation versus the WA75 and WA92 data. 



rric = 1.5 GeV, 

fip = 2/io, /ijj = /io, where 



/^o = J{pl + pl)/2 + ml 



(2.12) 



• the MRSA' [Martin95] proton parton densities, with the corresponding value of 

Af = 152 MeV, 

• the SMRS2 [Sutton92] set for the pion parton densities. 

Choosing smaller values for the charm mass, for example, one would get a broader 
theoretical curve. Alternatively, one could keep the default mass value, but choose 
/ifl = /io/2 for the renormalization scale, instead of our default choice ^^ = jiq. As 
shown in ref. [Frixione94] , this would lead to a broadening of the A0 distribution 
similar to the one caused by an average intrinsic transverse momentum with (/c^) = 
0.5 GeV2. 

The WA75 collaboration, and recently the WA92 collaboration, published in 
refs. [Aoki92a, Adamovich96a] the distribution of the transverse momentum of the 
heavy-quark pair. The theoretical prediction supplemented with a fc-p-kick with 
(fc^) = 1 GeV^ cannot reproduce the WA75 data, while it is in rough agreement 
with the WA92 measurement, as displayed in fig. |^. Unlike the azimuthal correlation, 
the p^iQQ) distribution is affected by fragmentation effects, since these can degrade 
the momenta of the quark and antiquark by different amounts. Fragmentation ef- 
fects also moderate the pair transverse momentum arising from gluon radiation or 
from an intrinsic parton transverse momentum. We have verified that at the end, at 
Ef, = 350 GeV, the fragmentation always tends to soften the p%{QQ) distribution. 
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Figure 8: NLO QCD result for tfie p%{QQ) supplemented witti an 
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Figure 9: Azimuthal correlation of DD pair versus the perturbative 
result in photoproduction for the E687 (left) and NA 14/2 (right) exper- 
iments. 



Summarizing, the available experimental results for azimuthal cc correlation in 
hadron-hadron collisions show a tendency to peak in the back-to-back region A0 = tt, 
but the peak is less pronounced than the one predicted by perturbative QCD. The 
addition of a kj^ kick of 1 GeV^ on average gives a satisfactory description of the data, 
if our default input parameters are chosen. The data on the p'^iQQ) distribution do 
not give a unique indication. While the theoretical prediction obtained with (/c^) = 1 
GeV^ is in rough agreement with the WA92 measurement, it is sizeably softer than 
the WA75 data. 
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Photoproduction of heavy quarks [Alvarez92, FrabettiQS, AdamovichST] is another 
example in which a kj, kick would induce broader A(j) and PriQQ) correlations. On the 
left-hand side of fig. P, the azimuthal correlation measured by the E687 collaboration 
is given, together with the NLO result. The NLO result supplemented by an intrinsic 
kr of the incoming partons is also shown, for {k"^) = 0.5 GeV^ and {k"^) = 1 GeV^. 
As can be seen, the data do not require a large intrinsic transverse momentum. All 
curves give a reasonable representation of the data, the one with (fc^) = 0.5 GeV^ 
being slightly better. A similar conclusion applies to the NA14/2 data (which are, 
however, affected by larger uncertainties), as shown on the right-hand side of fig. ^ 
The distribution in the transverse momentum of the heavy-quark pair is displayed in 
fig. |Ty. In this case, we see that the E687 data favour {k"^) = 1 GeV^. 
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Figure 10: Transverse momentum distribution of the DD pair versus 
the perturbative result for the E687 experiment. 



In conclusion, we have seen that the A0 and Pt{QQ) distributions are very sensi- 
tive to non-perturbative effects, especially in the hadropro duct ion case. We computed 
these distributions by assuming that the non-perturbative effects are parametrized by 
a fragmentation function and an intrinsic transverse momentum of the incoming par- 
tons. In general, the azimuthal correlation of the pair is well described by perturbative 
QCD supplemented by an acceptable /c^ kick. The recent result of ref . [Adamovich96a] 
on the Pt{QQ) distribution can also be described by the same choice of parameters. 
This is consistent with what was found in the study of the single-inclusive pr dis- 
tribution. However, the PriQQ) distribution reported by WA75 [Aoki92a] is much 
harder, and would require an unphysically large fcj, kick. 
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3 Heavy-flavour production at HERA 

The ep collider HERA offers new opportunities to study the production mecha- 
nism of heavy quarks and to test the predictions of perturbative QCD. The dominant 
contribution to the cross section is due to those events in which the virtuality of the 
photon exchanged between the electron and the proton is very small. In this case, the 
electron can be considered to be equivalent to a beam of on-shell photons, whose dis- 
tribution in energy (Weizsacker- Williams function [Weizsaecker34, Williams34]) can 
be calculated in QED. The underlying production mechanism is therefore a photo- 
production one, which has been studied extensively in fixed-target experiments. At 
HERA, the available centre-of-mass energy is about one order of magnitude larger 
than at fixed-target facilities (200 GeV versus 20 GeV). This energy regime is to- 
tally unexplored in photoproduction, and several new features have to be taken into 
proper account. In particular, the large contribution of the hadronic photon compo- 
nent introduces in the theoretical predictions a source of uncertainty that is normally 
negligible at fixed-target energies. 

A complementary way of studying heavy-flavour production at HERA is to retain 
only those events characterized by a large photon virtuality (DIS). Although the total 
rates are much smaller than the photoproduction ones, the hadronic component is 
completely eliminated, and more reliable theoretical and experimental results can be 
obtained. Also, the dependence of the data upon the photon virtuality can be used 
as a further test of QCD predictions. A significant improvement is expected with the 
planned luminosity upgrade of the HERA collider [Ingelman96] , when the number of 
charm (bottom) quarks produced will be of the order of 10^ (10^). 

In this section, we discuss some aspects of the perturbative QCD calculations of 
the heavy-flavour photoproduction cross sections of relevance for HERA. We review 
the phenomenology of charm and bottom production, comparing the available data 
with the theoretical predictions. We then discuss future perspectives. 

3.1 Photoproduction cross sections 

It is well known that an on-shell photon has a finite probability to fluctuate into a 
hadronic state before undergoing a hard collision. In this case, the photon is referred 
to as "hadronic" (or "resolved"), to contrast with those events in which it directly 
interacts with the hadron ("point-like" or "direct"). Therefore, a differential photon- 
hadron cross section can be written as the sum of a point-like and a hadronic photon 
contribution: 

rfa(^-)(P„ P,) = da^;j^,{P,, P,) + rfa[lJ)(P„ P,) . (3.1) 
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Thanks to the factorization theorems in QCD [Elhs79, Colhns89], we have 

j 

X daij{xP^, yPa, as(yu'^), /i'^, /i^, /i^) . (3.3) 



In eq. ( p.3|) , the /j are the partonic densities in the photon. Their physical meaning 
is analogous to the one of the more familiar partonic densities in the hadron, fj- 
They are universal, but not calculable in perturbation theory. They satisfy a renor- 
malization group equation [Witten77, Llewellyn 78, Bardeen79, DeWitt79] that can 
be obtained by slightly modifying the usual Altarelli-Parisi [Altarelli77] equation 

^-|^fl.+Si:i'««/r. (3.4) 

At the lowest order we have 

P,^ = N,e^^ [x^ + {l-xf) , (3.5) 

where A'^c = 3 is the number of colours and e^ is the electric charge of the parton i in 
units of the charge of the positron (for gluons, Cj = 0). The first term on the RHS 
of eq. (|3.4|) , which is not present in the evolution equation of the hadron densities, is 
due to the direct coupling of the photon to the quarks. 

To understand the strict interplay between the point-like and the hadronic compo- 
nent in eq. (|3.1|) , it is worth while to sketch the derivation of the subtracted partonic 
cross sections appearing in eqs. ( p.2| ) and (|3.3|). As a first step, one performs a direct 
calculation of the partonic cross sections for the processes •-} + j^Q + Q + X and 
i + j—yQ + Q + X (here i and j represent a generic parton). After the ultraviolet 
renormalization, one is left with the quantities da.yj and daij, which still display di- 
vergences due to the collinear emission of a massless final-state parton from one of 
the incoming partons. These divergences are subtracted by applying the prescription 
of the factorization theorem. Writing the expansion of the partonic cross sections at 
the next-to-leading order in QCD as 

da^j = aemC(sd^jj + aemalda!:^^ , (3.6) 

daij = a^dalj + alda^j , (3.7) 
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and an analogous expansion for the subtracted cross sections da^j and ddij, we get, 
using dimensional regularization (c? = 4 — 2e) and order by order in perturbation 
theory, 

da!^^{pi,P2) = da^^^{pi,p2), dalf{pi,p2) = dalf{pi,p2), (3.8) 

dcr!yj{pi,P2) = da^^]{pi,p2, -) + ^^ J dx i-Pk-yix) - Hk^{x)\ daf^{xpi,p2) 

+ ^J2jdx (ip,,(a;) -4f (a;)) da^^iip^,xp2), (3.9) 

k 

d^ff{pi,P2) = da^f{pi,p2, -) + ^^ J dx (-Pkiix) - kI^\x)] da^^j'{xpi,p2) 
+ ^T.Jdx {\PkA^)-Kif{x)^ dall\p„xp2). (3.10) 



On the RHS of eqs. (|3^ ) and (|3.1CI| ) we explicitly indicated the dependence of the par- 



tonic cross sections upon 1/e to recall that these quantities are still infrared- divergent. 
The divergences are nevertheless properly cancelled by the terms proportional to the 
Altarelli-Parisi kernels appearing in these equations. The functions K^" , K^^J and 
iffc7 are completely arbitrary, in that they define an extra finite part of the sub- 



traction; different choices correspond to different subtraction schemes. In eqs. (|37 
and ( |3.10| ) the MS scheme is equivalent to if = -ft' = 0. For greater generality, we 
have allowed the possibility to have different subtraction schemes on the photon and 
hadron legs. 

The key feature of eqs. (|3.9| ) and (|3.10|) is the term proportional to Pk-yix) in 



eq. ( ^.91 ). The physical origin of this contribution is the direct coupling of the photon 
with the quarks, which in turn implies the presence of the inhomogeneous term in 
eq. (p.4|) . In the language of the Altarelli-Parisi equations, this means that the 
infrared divergences due to the collinear emission of quarks from the incoming photon 
in the point-like component, eq. ( |3.2| ), are re-absorbed into the partonic densities in 
the photon, which appear in the hadronic component, eq. (|3.3 )- The same argument 
can be formulated in terms of renormalization group (RG) equations. Equation (|3.2| ) 
is RG- invariant with respect to the variation of the scales /x^ and fip] eq. ( |3.3| ) is RG- 
invariant with respect to the variation of the scales fi'^ and fi'p. But neither eq. (|3.2| ) 
nor eq. (|3.3|) are separately RG-invariant with respect to the variation of the scale /x^; 
when varying /i^ in eq. (P73|), a residual dependence is left, due to the inhomogeneous 
term in eq. (|3^ ), which is cancelled only by the explicit dependence of eq. (|3.2|) upon 
/x^. A third way of understanding this issue is to consider a change of the subtraction 
scheme. The partonic densities are modified as follows: 

ir' = /f^+|E<^®/r^ (3.11) 

i 

j 



Heavy-flavour production at HERA 23 



The term Hi^, which is due to the change of scheme of the photon densities entering 
eq. (|3.3| ), affects eq. ( |3.2|) , since it defines the finite part of the infrared subtraction 
in the photon leg via eq. ( W^ - 

It should now be clear that the point-like and the hadronic components of the 
photoproduction cross sections are very closely related, and that only their sum is 
physically meaningful. The separation of a cross section into a point-like and a 
hadronic component is ambiguous beyond leading order (for a detailed discussion, 
see refs. [Schuler93, Schuler93a]), as is shown by the discussion on the change of 
subtraction scheme: finite terms can be subtracted to one piece and added to the 
other, without affecting physical predictions. One explicit example is presented in 
ref. [Gluck92], where a factorization scheme (DIS^) for the photon densities is in- 
troduced, which uses K^'^^ = and H ^ 0. That notwithstanding, we will keep 
on talking about the point-like and the hadronic component. The reason for this is 
twofold: we can use a leading-order approximation to get a physical picture of the 
separation between the two components; and the term (H) that can be exchanged 
between the two components is numerically small. 

The photon parton-densities are quite soft. Therefore, the hadronic component 
is only important for large CM energies and small masses for the produced system. 
We will see in the following that it potentially affects charm production at HERA. 
This process can therefore be used to constrain the densities in the photon, which are 
experimentally very poorly known at present. 

As a final step, in order to obtain the QQ cross sections in electron-proton col- 
lisions, the photoproduction cross sections must be convoluted with the Weizsacker- 
Williams distribution: 






f^'\y) = ^ ^^ ^"^ iog ^-V7' • (3-13) 



The scale fiw^ is, in general, a function of y determined by the kinematics of the 
process considered. It has been pointed out in refs. [Budnev74, 01sen79, Bawa89, 
Catani91, Frixione93] that the choice of fi^w should also take into account the dy- 
namics of the production mechanism. In the case of the production of heavy quarks it 
is reasonable to set ^ww = ^"^q, where ^ is of order 1. It was shown in ref. [Frixione93] 
that for the point-like component an appropriate choice for the parameter ^ is .^ = 1. 
The hadronic component case is more involved, since the partonic densities in the vir- 
tual photon introduce in the problem an additional mass scale. The presence of this 
scale is due to the fact that the densities rapidly fall to zero when the virtuality of the 
photon approaches the hard scale of the process [Uematsu82, Borzumati93, Drees94], 
which is of the order of the heavy-quark mass. It follows that the scale fiww niust 
be chosen smaller than mg. In ref. [Frixione95a] it has been argued that ^ = 0.6-0.7 
gives sensible results. Although this is only a rough estimate, the uncertainty due to 
the choice of the parameter ^ is much smaller than the uncertainties coming from the 
phenomenological parameters entering the calculations. 
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It is interesting to notice that a photoproduction event can be experimentally 
defined by means of an anti-tag condition: all those events in which the electron 
is scattered at an angle larger than a given value 6c are rejected. The form of the 
Weizsacker-Williams function in the presence of an anti-tag condition can be explic- 
itly worked out [Frixione93] . Indicating with E the energy of the incoming electron 
in the laboratory frame, we get 



/r(i/) 




mly"^ 



(3.14) 



The error due to the approximation can be estimated to be 0{6l, ml/E"^) in this case, 
and it is therefore quite small for applications to HERA physics, where 6'c ~ 5 x 10~^. 
Also notice that the non-logarithmic term is singular in y and therefore represents a 
non-negligible correction. 

3.2 Charm photoproduction 

Taking into account all the relevant sources of uncertainty, which have already 
been discussed in section ^, we get a theoretical prediction for the total rates with 
a large error. At the NLO, the point-like cross section changes by a factor of about 
4 when varying the mass in the range 1.2 GeV < nic < 1.8 GeV, and by a fac- 
tor of 2 when varying the renormalization scale. The choice of the proton partonic 
densities induces a 50% uncertainty at JS^p = 30 GeV, and a factor of about 5 

[Frixione95b] at JS^p = 300 GeV. These effects are even more dramatic in the case 
of the hadronic component, where the dominant source of uncertainty is the choice of 
the partonic densities in the photon. Results from dijet photoproduction at HERA 
give some indication on the gluon density in the photon [Ahmed95], but the very lim- 
ited statistics does not allow either a distinction between different NLO parametriza- 
tions [Gluck92a, Gordon92, Aurenche92, Aurenche94, Gordon96], or a discrimination 
between the LACl [Abramowicz91] gluon and the flatter ones of the NLO sets. Other 
LO sets [Hagiwara95] are also consistent with the data. A reasonable way to estimate 
the effect due to the uncertainty in the gluon density of the photon is to take the sets 
GRV-HO [Gluck92a] and LACl as the two extremes. At the highest photon-proton 
centre-of-mass energies available at HERA, the prediction for the hadronic compo- 
nent with LACl is one order of magnitude larger than the prediction with GRV-HO, 
and much larger than the point-like component as well. 

The comparison between the theoretical predictions and the experimental results 
[Derrick95, Aid96] is presented in fig. ^, where only the uncertainty on the QCD 
result due to the choice of the renormalization scale is displayed. Despite the fact 
that, as discussed before, this only partially accounts for the overall uncertainty on 
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Figure 11: Total cross section for tfie ptiotoproduction of cc pairs, as 
a function of ttie 7p centre- of -mass energy: next-to-leading order QCD 
predictions versus experimental results. 



the theoretical predictions, the data show a satisfactory agreement with the curves 
obtained by choosing nic = 1.5 GeV and MRSG and GRV-HO for the partonic densi- 
ties in the proton and the photon respectively. It is important to notice that a single 
choice of the input parameters allows us to describe the data in the whole energy 
range considered; this is a non-trivial result, the energies available at HERA being 
one order of magnitude larger that those available at fixed-target experiments. The 
importance of the hadronic component of the photon at fixed target and at HERA is 
also clear from fig. |TT|. 

The experimental determination of the total charm cross section at HERA deserves 
some further comment. The experiments are sensitive to production in the central 
rapidity region, typically \y\ < 1.5, where the cross section is far from its maximum 
(see fig. ^). Furthermore, a small-p^ cut is applied to the data, in order to clearly 
separate the signal from the background. In order to get the total cross section, one 
has therefore to extrapolate to the full rapidity and transverse momentum range. 
Large rapidities and small transverse momenta typically involve small-x values, and 
the extrapolation is therefore subject to uncertainties due to our ignorance of the 
small-x behaviour of the parton densities. Since these uncertainties are very large, 
we believe that it would be useful to present also the measurement of the cross 
section limited to the directly accessible rapidity region. The lack of sensitivity to 
large negative rapidities is one of the reasons why it is at present impossible to 
use the photoproduction of charm quarks at HERA to distinguish among different 
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Figure 12: Transverse-momentum distribution of D* mesons in ep 
collisions at HERA, for Q^ < 0.01 GeV^ . Experimental results [Aid96] 
are compared with next-to-leading order theoretical predictions. 



parametrizations of the gluon density in the proton. An equally important reason (as 
pointed out in section H) lies in the presence of inconsistencies in the low-energy data. 
From fig. |ll] no definite conclusion can be drawn on the photon densities either. 

The first measurement of single-inclusive distributions for charm photoproduction 
at HERA has been presented by the HI Collaboration [Aid96]. The data are inte- 
grated over a large range of photon-proton energies, and the results should therefore 
be compared with the theoretical predictions for electroproduction in the Weizsacker- 
Williams approximation. 

Since the Weizsacker-Williams function grows in the small-x region, the bulk 
of the contribution to the cross section is due to photons of relatively low energy. 
Therefore, the role of the hadronic component of the photon and of the small-x 
region in the proton densities is marginal. It turns out that, when cuts similar to 
the experimental ones are applied and electroproduction in the Weizsacker-Williams 
approximation is considered, the shape of the theoretically predicted distributions is 
very stable [Frixione95a] . It follows that single- inclusive charm electroproduction is 
of little help in constraining the partonic densities of both the proton and the photon, 
but can be used as a valuable test of the production mechanism. 



Figure |T2| presents the comparison between theory and HI data [Aid96] for the p^ 
spectrum of the D* mesons. The data are relative to the so-called "tagged" sample, 
defined by the condition Q^ < 10~^ GeV^. The QCD predictions for charm quarks 
have been convoluted with the Peterson fragmentation function [Peterson83]. Given 
the limited statistics of the measurement, the comparison of the data with the theory 
appears to be satisfactory, and qualitatively analogous to similar results at fixed- 
target experiments [Frixione94] . A comparison has also been carried out by the HI 
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Figure 13: Left: rapidity distribution of D* mesons in ep collisions at 
HERA, for Q"^ < 0.01 GeV^ . Experimental results [Aid96] are compared 
with next-to -leading order QCD predictions. Right: pseudorapidity dis- 
tribution of fragmented charm quarks in monochromatic photon-proton 
collisions. 



Collaboration [Aid96] for the untagged sample ((Q^) — 0.2 GeV^); the agreement 
is slightly worse, but still satisfactory. As shown in fig. [T^, the uncertainty on the 
theoretical prediction due to the value of the charm mass is completely negligible 
when compared to the uncertainties in the experimental data. 

The comparison between theory and HI data for the rapidity distribution is pre- 
sented on the left-hand side of fig. 0. Although the agreement is acceptable, only 
one point being more than one standard deviation away from the theoretical predic- 
tions, higher statistics are required in order to perform a more significant comparison. 
Indeed, preliminary HI results^ show that, when improving the statistics, the data 
point at y{D*) = 0.25 is also close to the theoretical curve. We point out, however, 
that recent preliminary data from ZEUS [Zeus96a] are not in satisfactory agreement 
with QCD expectations. 

If data with larger statistics were available, it would be possible to consider the 
production processes initiated by very energetic photons only. The right-hand side 
of fig. [1^ shows the theoretical prediction for the pseudorapidity distribution of frag- 
mented charm quarks in monochromatic photon-proton collisions, at a centre-of-mass 



energy of 



286 GeV, and for different choices of the partonic densities in the 



photon. As is apparent from the figure, the results obtained with the LACl and 
GRV-HO sets are completely different in the forward region. Even at the moderate 
pseudorapidity values covered by the present configuration of the detectors, the large 
difference induced by the two photon sets should have measurable effects. 



^C. Grab, private communication. 
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Because of the large photon-proton centre-of-mass energy available at HERA, the 
log(5'/m^) terms appearing in the charm cross section may get large and spoil the con- 
vergence of the perturbative series. The problem of resumming these terms {small-x 
effects) has been tackled by several authors [EllisQO, Catani89, Catani90, CollinsQl], 
mainly in the context of b production in hadronic collisions. In ref. [Frixione95b] an 
estimate has been given of the importance of the small-x effect for charm physics in 
the HERA energy range. It has been found that, resumming the log(S'/m^) terms, 
the total rate can be increased by by 20% to 40% with respect to the next-to-leading 
order prediction for the point-like contribution, and by 20% to 45% for the hadronic 
contribution, depending on the partonic densities used in the calculation (as a rule 
of thumb, the flatter is the gluon density at small x, the larger is the contribution 
expected from the resummation of the small-x effects). These effects are therefore 
negligible with respect to the other sources of uncertainty on the theoretical predic- 
tions. 

We finally observe that the transverse-momentum distribution is in principle af- 
fected by the presence of log{pT/mc) terms. These logarithms can be resummed by 
observing that, at high Pt, the heavy-quark mass is negligible, and by using per- 
turbative fragmentation functions [Cacciari96, Kniehl95]. Remarkably enough, the 
fixed-order and the resummed results of [Cacciari96] agree in a very wide range in 
Pt- Recently, the massless approach has been used in ref. [Kniehl96] to predict the 
D* spectrum at HERA. 

3.3 Bottom photoproduction 

Thanks to the higher value of the quark mass, perturbative QCD predictions for 
bottom production are more reliable than those for charm. In fact, all the uncer- 
tainties we have discussed for charm are in this case strongly reduced. The resum- 
mation of the small-x effects has been estimated in ref. [Frixione95b] to increase the 
next-to-leading order result by ~ 5%. In monochromatic photon-proton collisions, 
the point-like component has an uncertainty of a factor of 2 if all the parameters 
are varied together in the direction that makes the cross section larger or smaller. At 
JS-yp = 100 GeV, the lower and upper limits of the point-like component are 16 nb and 

35 nb respectively, while at JS^p = 280 GeV we get 41 nb and 101 nb [Frixione95b] . 
The hadronic component has larger uncertainties, but much smaller than in charm 
production. As discussed above, the main source of uncertainty is the parton density 
set for the photon. Nevertheless, in bottom production the small-x region is probed 
to a lesser extent than in charm production, and the sensitivity of the result to the 
photon densities is therefore milder; we get an uncertainty of a factor of 3 (to be 
compared with a factor of 10 in the case of charm). The hadronic component can 
still be the dominant contribution to the photoproduction cross section, if the gluon 
in the photon is as soft as the LACl parametrization suggests. 
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Figure 14: Full uncertainty on tfie transverse-momentum distribution 
for bottom electroproduction (Weizsdcker- Williams approximation) with, 
Peterson fragmentation and a pseudorapidity cut. 



The bottom rates are about a factor of 200 smaller than the charm ones. To 
perform a statistically significant study of bottom production, the luminosity up- 
grade at HERA is necessary. In any case, it is very likely that a comparison with 
the theory could only be done by considering electroproduction in the Weizsacker- 
Williams limit. In this case, consistently with what was discussed for charm, the 
sensitivity of the theoretical predictions to the input parameters is sizeably reduced, 
and a reliable comparison between theory and data can be performed. For example, 
in electroproduction the hadronic component contribution to the total cross section 
is at most 75% of the point-like contribution, even if the LACl set is used. The most 
interesting results are, however, obtained when considering more exclusive quantities, 
such as single-inclusive distributions. In particular, as shown in fig. 0, the transverse 
momentum of the bottom quark at HERA can be predicted by perturbative QCD 
quite accurately. It is clear that even with the LACl set the hadronic component 
affects the prediction only marginally; this fact is basically a consequence of the ap- 
plied pseudorapidity cut. We can therefore regard fig. |14| as a reliable prediction of 
QCD for the Pt spectrum of h hadrons at HERA. The comparison of this prediction 
with the data would be extremely useful in the light of the status of the comparison 
between theory and data for h production at the Tevatron (see section Hf). 



3.4 Deep-inelastic production 

The first data on charm production in the deep-inelastic regime at HERA have 
recently become available [Zeus96, Adloff96]. On the theoretical side, NLO QCD cal- 
culations for total and single-inclusive deep-inelastic production have been performed 
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Figure 15: Comparison between tfieoretical predictions and experimen- 
tal data for single-inclusive distributions in DIS charm production. 



in refs. [Laenen92a, Laenen93, Laenen93a]. Recently, fully exclusive cross sections 
have been computed [HarrisQS, Harris95a]. In this context, it is possible to perform 
comparisons between theory and experiment in different aspects. 

To start with, one can consider the contribution of charm quarks to the proton 
structure functions F2 and F^. It turns out that NLO corrections to these quantities 
are non- negligible. The full NLO results display a mild dependence upon the renor- 
malization and factorization scales (well below 10% for Q"^ larger than 10 GeV^, when 

varying scales between half and twice the default value 



/Q^ + m"^). The dominant 
uncertainty in the calculation is the small-x behaviour of the gluon density in the 
proton. For this reason, it is possible to conclude that structure-function measure- 
ments may provide useful information about the gluon distribution in the proton (see 
ref. [Laenen96a] for a discussion of this issue). For Q"^ ^ 4m^ the calculation of 
F2 and Fl requires the inclusion of large logarithmic effects arising from evolution. 
Approaches to this problem can be found in refs. [Aivazis94, Aivazis94a, Lai97] and 
[Martin96]. 

In fig. |I5|, taken from ref. [Laenen96a], we show experimental data for some single- 
inclusive distributions, superimposed on QCD theoretical predictions at NLO. The 
following distributions are shown: a) the D* transverse momentum p!^* , b) its pseu- 
dorapidity rj^' , c) the hadronic final state invariant mass and d) the Q^ one, for the 
kinematic range 5 GeV^ < Q^ < 100 GeV^, 1.3 GeV < p^* < 9 GeV and \r]^*\ < 1.5. 
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A cut on the DIS y variable < y < 0.7 is also applied. The D* momenta are ob- 
tained by applying a Peterson fragmentation function to the produced charm quark. 
Also shown are experimental data of the ZEUS collaboration [Zeus96]. 

The shape of the distributions considered is significantly modified by the intro- 
duction of radiative corrections. As we can see, the data are in reasonable agreement 
with theoretical predictions. More statistics is needed for a meaningful comparison 
with the theory. 

3.5 Future physics 

In this subsection we will discuss a few topics on future physics possibilities in 
heavy-fiavour production at HERA. With the planned upgrades of the HERA collider, 
an integrated luminosity of 100 pb~^ or larger will be achieved. This will probably 
allow the production of a sizeable sample of double-tagged charm events [Eichler96] , 
and therefore to study double-differential cross sections. The possibility of having 
polarized beams has also been considered [Ingelman96] . Furthermore, the HERA-B 
fixed-target program will become operational. 

3. 5. 1 Determination of f^^' 

An interesting application of double-differential measurements is the determination 
of the proton gluon density /j^^ [BuchmuellerQl, Riemersma92, Frixione93a] . At 
present, this quantity is not directly measured: DIS data allow its indirect extraction 
from the Q^ evolution of the structure functions, and direct photon and inclusive jet 
measurements constrain it in complementary regions of x and Q^. 

In photon or electron-hadron collisions, the advantage is that the gluon density 
enters the cross section in a simpler way than in the case of hadronic collisions. It is 
easy to show that this in turn implies that, if one could determine the invariant mass 
and rapidity of the produced system, the cross section would be directly proportional 
to flj^^ and to other calculable factors. For example, the LO heavy-lavour electropro- 
duction cross section can be written, in the Weizsacker-Williams approximation, and 
neglecting the hadronic component, in the following way: 

"^""^'^ ^^J^ = 4 /f (^.>^o)/f (^.,/^l)^SnM^«), (3.15) 



where Mqq is the invariant mass of the heavy-quark pair, and yq-Q is the rapidity of 
the pair in the electron-proton centre-of-mass frame (we choose positive rapidities in 
the direction of the incoming proton); E = ^/S is the electron-proton centre-of-mass 
energy, and 

x^ = -|^exp(-|/QQ), (3.16) 

Xg = -^expiygg). (3.17) 
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Figure 16: Xg distribution in ep collisions (Weizsdcker-Willianis ap- 
proximation) at HERA, for rric = 1.5 GeV and ttie MRSG proton par- 
ton- densities. 



The function /^^^^ is the Weizsacker-Wilhanis function we already discussed in sub- 
section p?T[ We observe that all quantities on the right-hand side of eq. ( |3.15| ) are 
calculable, except for /^^\ which can therefore be measured. The inclusion of ra- 
diative corrections does not pose any problem. This issue is discussed in detail in 
ref. [Frixione93a] , where it is shown that f^^^ can be measured up to NLO accuracy 
by using double-differential charm data. 

As discussed previously, the heavy-flavour cross section receives a large contribu- 
tion from the hadronic component. In order to extract the gluon density of the proton 
with the method outlined above, we have to consider only those kinematical regions 
where the hadronic component is suppressed. We study this possibility in flg. |16 , 
where we present the NLO QCD distribution in the variable Xg, in electron-proton 
collisions at yS = 314 GeV. The partonic densities in the proton are given by the 
MRSG set, while we considered both the LACl and GRV-HO sets for the photon. 



Figures |T6|a)-|T6|c) show the effect of applying a cut on the invariant mass of the pair. 
Even in the case of the smallest invariant-mass cut, there is a region of small Xg where 
the hadronic component is negligible with respect to the point-like one. When we in- 
crease the invariant-mass cut, we notice that the hadronic component decreases faster 
than the point-like one. This is due to the fact that, for large invariant masses, the 
production process of the hadronic component is suppressed by the small value of the 
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gluon density of the photon at large x. By pushing the invariant-mass cut to 20 GeV, 
it turns out that the point-hke component is dominant over the hadronic one for Xg 
values as large as 10~^. We then conclude that the theoretical uncertainties affecting 
the charm cross section, in the range of 10~^ < Xg < 10~^, are small enough to allow 
for a determination of the gluon density in the proton by using invariant-mass cuts 
to suppress the hadronic component. In a more realistic configuration, as the present 



one of the detectors at HERA, additional cuts are applied to the data. Figure |T6|d) 
shows the effect, on the Xg distribution, of a small-pr and a pseudorapidity cut ap- 
plied to both the charm and the anticharm. In this case, even without applying an 
invariant-mass cut, the point-like component is dominant in the whole kinematically 
accessible range. 

3.5.2 Polarization asymmetries 

It is conceivable that, in the future, the HERA collider will be operated in a po- 
larized mode. At present, a good result for the polarization of the positron beam 
has been obtained, and feasibility studies for the polarization of the proton beam 
are under way [Ingelman96] . At leading order in QCD, the heavy-flavour production 
cross section in polarized ep collisions is proportional to the polarized gluon density 
in the proton, Ag. Therefore, data on charm production could be used to measure 
Ag directly, as previously shown for the unpolarized case. In practice, the situation 
for the polarized scattering is much more complicated. First of all, it is not possi- 
ble to perform a complete NLO analysis, because NLO cross sections for polarized 
charm photoproduction have never been computed. Furthermore, there is no exper- 
imental information on parton densities in the polarized photon. It is reasonable, 
however, that charm production at the HERA collider in the polarized mode can 
help in constraining the polarized gluon density in the proton. This possibility was 
first suggested in refs. [Gluck88, GluckQl, VogelsangQl], and recently reconsidered in 
refs. [Frixione96b, Stratmann96] . 

In order to reduce the impact of the uncertainties induced by radiative corrections, 
it is useful to present predictions [Frixione96b] for the ratio Ao jo (asymmetry), where 
a is the unpolarized cross section and 



Ag 



l(att-aTi). (3.18) 



Here a^^ and a^^ are the cross sections for cc production with parallel and antiparallel 
polarizations of the incoming particles respectively. One might expect that the effect 
of the radiative corrections approximately cancels in the ratio. It must be stressed 
that, for consistency, the unpolarized cross section o appearing in the asymmetry 
must be calculated at the leading order, as the polarized one. The NLO value of a 
can then be used to estimate the sensitivity of the experiment. 
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Figure 17: Asymmetry cross section versus transverse momen- 

tum in polarized ep collisions (Weizsdcker- Williams approximation) at 
\/s = 314 GeV. The minimum observable asymmetry, computed at 
next-to-leading order, is also displayed. 



A rough estimate of the minimum value of the asymmetry observable at HERA 
can be obtained by requiring the difference between the numbers of events with par- 
allel and antiparallel polarizations of the initial-state particles to be larger than the 
statistical error on the total number of observed events. This gives 



Aa 



a 



1 



V2(T£e' 



(3.19) 



where C is the integrated luminosity and the factor e accounts for the experimental 
efficiency for charm identification and for the fact that the colliding beams are not 
completely polarized. This procedure can be applied to total cross sections, as well as 
to differential distributions; in this case, the values of a and Aa have to be interpreted 
as cross sections per bin in the relevant kinematical variable. 

In ref. [Frixione96b] it was shown that total cross section asymmetries for the 
point-like component are quite small in absolute value, and can be measured only if e 
is equal to or larger than 1% (0.1%), assuming C = 100 pb~^ (1000 pb~^). Therefore, 
even with a vertex detector (see ref. [Eichler96]), it appears unlikely that this kind 
of measurements will be performed at HERA. Furthermore, in ref. [Stratmann96] a 
rough estimate of the hadronic contribution to the polarized cross section has been 
given, assuming polarized parton densities in the photon to vary between zero and the 
corresponding unpolarized densities in order to get a lower and an upper bound on 
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the cross section. It was found that a non-neghgible contamination of the point-hke 
resuh might indeed come from the hadronic process. The situation clearly improves 
when considering more exclusive quantities; in ref. [Frixione96b] it was found that at 
moderate pr values the asymmetry for the point-like component can be rather large, 
well above the minimum observable value (in this region, the experimental efficiency 
is sizeable [Eichler96]); this is shown in fig. [l^. In ref. [Stratmann96] it was argued 
that the hadronic component should have a negligible impact in this case. 

We conclude that charm data in high-energy polarized ep collisions will help in 
the determination of the polarized gluon density of the proton, provided that the 
integrated luminosity and the experimental efficiencies will allow a study of at least 
single-inclusive distributions. 

3.5.3 HERA-B 

The HERA-B program will allow a detailed study of h production in proton-copper 
collisions, at v^ = 39.2 GeV. Since this energy is relatively close to threshold, one 
may worry about the presence of large soft gluon resummation effects. Large resum- 
mation effects were found in ref. [Kidonakis95]. In ref. [Catani96a] it was instead 
argued that these large effects have a spurious origin, and a properly performed re- 
summation gives corrections of the order of 10%. This issue will be discussed in 
more detail in section ^ For example, from the NLO calculation with the MRSA' 
[Martin95] parton densities and nib = 4.75 GeV, we get a^b = 10.5 ^^'.j, nb, a range 
obtained by varying the renormalization and factorization scales from ?Tib/2 to 2mi,. 
Thus the upper band is 80% higher than the central value, to be compared with a 10% 
increase from the resummation effects of ref. [Catani96a]. This result is much less 
dramatic than the results of ref. [Kidonakis95] . For completeness we report in table ^ 
cross section values obtained using different sets of structure functions, corresponding 
to different values of A^^, and to different values of the bottom quark mass. Observe 
that measurements of h cross sections in pN collisions at a beam energy of 800 GeV, 
close to the HERA-B regime, have been presented in refs. [Jansen95, Spiegel96]. 
These measurements (see fig. in section ^ are consistent with the range reported 
in table |[ 

4 Heavy-quark production at hadron colliders 

Hadron collisions are by far the largest source of heavy quarks available today. 
While the environment of high-energy hadronic interactions does not allow to trig- 
ger on the largest fraction of charm and bottom produced, the production rates are 
so large that the number of recorded events allows today 6-physics studies that are 
competitive with those of e^e~ experiments. The introduction of new experimen- 
tal techniques, such as the use of silicon vertex detectors, which enable the tagging 
of events containing bottom quarks, led in the recent years to high-statistics and 
low-background measurements of the 6-production properties over a large range of 
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Table 2: Total cross sections in nb for bottom production at HERA-B. 
The central values correspond to fi = nib, upper and lower values to 
H = mb/2 and fi = 2nib respectively. The structure functions are those 
of refs. [Martin95, Martin95a]. 



transverse momenta. The high energies and luminosities available at the Tevatron pp 
collider, in addition to the enhanced 6-tagging capability, have furthermore led to the 
discovery of the top quark [Abe94, Abe95, Abachi95], opening the way to a new set 
of tests of the heavy-quark production dynamics. 

From the point of view of QCD studies, heavy-flavour production in high-energy 
hadronic collisions has better potentials than in fixed-target experiments. The b 
quarks produced at large pj. can be studied in perturbative QCD with smaller contam- 
ination from non-perturbative effects. For example, the impact of the initial-parton 
transverse momentum is largely reduced with respect to fixed-target charm produc- 
tion. Furthermore, the fragmentation properties of heavy fiavours in high transverse 
momentum jets can be directly studied, since the transverse momenta involved are 
typically perturbative. Finally, top production should be fully perturbative, and 
therefore the ultimate testing ground for the theory of heavy-fiavour hadroproduc- 
tion. 

In this section we review the current status of these measurements and their 
implications for the study of QCD. We discuss here first the phenomenology of charm 
and bottom production. We cover the inclusive production of b quarks, b hadrons, 
jets containing b and c quarks, bb correlations, and we finally discuss the associate 
production of heavy fiavours and electroweak gauge bosons. We then discuss the 
phenomenology of top-quark production. The large mass of the top quark implies that 
theoretical predictions based on NLO calculations should provide a rather accurate 
description of its inclusive production properties, with uncertainties in the range of 
10%. The experimental statistics currently available cannot probe the theory at this 
level of accuracy, and significant improvements will only occur with the next round 
of data taking at the Tevatron and with the LHC. Nevertheless, the current data 
are an important step forward, and it is encouraging that they basically confirm the 
theoretical expectations. 
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4.1 Inclusive bottom production 

4-1.1 Preliminaries 

The distribution most commonly studied by the hadron-collider experiments is the 
6-quark differential pr spectrum, integrated within a fixed rapidity range and above 
a given p^ threshold (p"'" ) : 

72 

a(pr)= / dy f . dpr-^. (4.1) 

J\y\<yma^ JPt>pT" dydpr 

The UAl experiment at the CERN pp collider used ymax = 1-5, while the CDF and 
DO experiments at the Tevatron use ymax = 1- QCD calculations for this distribution 
are available at NLO [Nason89, BeenakkerQl]. 

In analogy to the case of fixed-target charm production, the study of NLO b cross 
sections indicates a strong dependence on the choice of factorization and renormal- 
ization scales. We show in fig. [1^ this dependence at various values of Pt for the 
two energies y/S = 630 and 1800 GeV. Notice the growth of the cross sections at 
small scale values. Notice furthermore that for values of p"'" < 40 GeV the scale 
dependence at 630 GeV and 1800 GeV are approximately the same. The large scale 
dependence indicates that corrections of higher orders are significant both at 630 and 
at 1800 GeV. Above p"'" = 40 GeV, on the other hand, the scale dependence at 
630 GeV is significantly reduced, suggesting a more reliable perturbative expansion. 

In the following we will present results obtained by varying simultaneously /xr and 
/ip within the range0/io/2 < /i < 2/io (/ig = Pt'^ + ^), which represents an acceptable 
window within which to explore the scale dependence. 

The uncertainty induced by the choice of different sets of parton density functions 
(PDF) is rather marginal, given the tight constraints set by the DIS and HERA data. 
One should, however, remember that the gluon density, which is the most important 
ingredient in b and c production at colliders, has not yet been measured directly 
for the values of x probed by the b measurements at CERN and at the Tevatron, 
approximately limited to the range 0.01 < x < 0.1. As a standard set of parton 
densities we will choose here the MRSA' set [Martin95], for which Af^ = 152 MeV. It 
should be pointed out, however, that this value of A|^^ yields a value of ag significantly 
lower than that extracted from different observables [Barnett96] . Parton distribution 
fits have therefore also been performed with fixed values of A|*^ [Martin95a, Lai96]. To 
explore the dependence of our results on the choice of as, we will also present results 
obtained using the PDF set MRS125 [Martin95a], which was extracted by forcing 
as{Mz) = 0.125. Notice that this procedure is not necessarily consistent, since the 
inconsistency in the extraction of ag is present in the data themselves. Should the 
future DIS measurements of ag become closer to the LEP value, as suggested by recent 
results [Harris96, Seligman97, Schmelling96], it is very likely that also the PDF fits, 
and with them our results, will change. 

Given the large values of pr probed by the collider experiments, the 6-mass de- 
pendence of the theoretical result is small. We will choose as a default the value 
rrib = 4.75 GeV, and allow for variations in the range 4.5 GeV < mi, < 5 GeV. 

^Unless otherwise specified we will always assume the renormalization and factorization scales to 
be equal, and will simply denote them by fi. 



38 Heavy-Quark Production 



6^ 



a 

b 



3.0 




Figure 
at 630 



18: Scale dependence of the integrated b-quark p^ distribution 
GeV (dashed lines) and at 1800 GeV (solid lines), for different 



values ofp^^". 



4-1.2 The effect of higher- order corrections 

The significant scale dependence of tlie NLO results for bottom production is a symp- 
tom of large contributions from higher orders. We briefly discuss here two possible 
sources of higher-order corrections: small-x effects, possibly relevant in the low-pj, 
domain, and large logarithms, which appear at high Pt- 

The possibility that higher-order corrections at the Tevatron may be very large 
had been pointed out early on in ref. [Nason88]. One can show that, in the high- 
energy limit (p = Amq/S -^ 0), large corrections arise to all orders in perturbation 
theory. The source of these large corrections is the presence, starting from 0{af), of 
diagrams where a gluon exchanged in the t channel becomes soft. In this limit, it is 
easy to show that the dominant behaviour of the Born and of the 0{a^) term are 



o-Born{gg^QQ) OC -^ p , 



(4.2) 



^i3){99^QQ) 



OC 



a 



m 



(4.3) 



where p = Arn^/ s. A crude but revealing estimate of the relative size of leading- 
order and higher-order terms can be obtained by assuming a simple form of the gluon 
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densities at small x, e.g. /(x) = Ajx^^^ , with 5 < 1. It is easy to show that the total 
cross sections scale as follows: 

n, . (nn^nn\ "^ log - (if 5 log 1/p < 1) 

asornigg^QQ) [a,i±^ (if <5logi/p > i) 

Notice that in the Feynman scaling approximation (e.g. x g{x) constant at small x), 
the 0{al) correction is enhanced by a large logarithm of the ratio between the total 
hadronic CM energy and the heavy-quark mass. Logarithms of this type, also known 
as small-x logarithms, arise at all orders of perturbation theory, and need to be re- 
summed. The resummation theory for this class of corrections has been applied to the 
problem of heavy-flavour production in refs. [Ellis90, CollinsQl, CataniQO, Catani91, 
Levin91]. The results contained in ref. [CollinsQl] indicate an enhancement of the 
NLO total bottom cross section at the Tevatron by no more than 30% due to the 
resummation of small-x effects, the highest values being obtained with small-5 gluon 
densities. Naive reasoning based upon eq. ( |4.4| ) would lead to a similar pattern. 



These results suggest that small-x effects at the current collider energies are not 
large enough to justify discrepancies as large as those initially found at the Teva- 
tron. Use of a gluon density more singular than 1/x, such as for example recent 
parametrizations obtained from the HERA data, should reduce these effects even 
more. 

Another class of potentially large higher-order corrections appears when h quarks 
are produced at high transverse momentum. At large momentum, in fact, the h quark 
behaves more and more like a massless particle, radiating an increasingly large amount 
of its energy in the form of hard, collinear gluons. This physical phenomenon is 
associated with the presence of logarithms of the form log(pT/m), which appear at all 
orders in perturbation theory. This problem was already examined in ref. [Nason89]. 
Some logarithmically enhanced higher-order terms were estimated, and found to be 
negative and small up to transverse momenta of the order of 20 GeV. 

Techniques are available to resum this class of logarithms. Two different groups 
have approached this problem in the recent past. Cacciari and Greco [Cacciari94] 
have folded the NLO cross section for production of a massless parton i {i = g,q) 
[Aversa89] with the NLO fragmentation function for the transition i^b [Mele90, 
Mele91]. The evolution of the fragmentation functions resums all terms of order 
a'2 log"' {pt / m) and a'g~^^log"'{pT/m). All the dependence on the 6-quark mass lies 
in the boundary conditions for the fragmentation functions. This approach ensures 
a full NLO accuracy, up to corrections of order m?/{m? + p\). In particular, this 
formalism describes NLO corrections to the gluon splitting process, which in the 
0{a1) calculation is only included at the leading order. One can verify, by looking 
at the Born- level production process as a function of the quark mass, that, in order 
for the mass corrections not to exceed the 10% level, it is necessary to limit the 
applications of this formalism to the region of p^ > 20 GeV. 



40 Heavy-Quark Production 



Figure |I9| shows the differential 6-quark p^ distribution obtained in the fragmentation- 
function approach, compared to the standard fixed-order NLO result. Several features 
of this figure should be noticed. To start with, the scale dependence is significantly re- 
duced with respect to the fixed-order calculation. Furthermore, in the range of appli- 
cability of this calculation (i.e. p^. > 20 GeV) the result of the fragmentation-function 
approach lies on the upper side of the fixed-order NLO calculation. The resummed 
calculation is, however, always within the uncertainty band of the fixed-order one. 
Finally, notice that the overall effect of the inclusion of higher-order logarithms is a 
steepening of the pj, spectrum, as is natural to expect. 

Another calculation has recently appeared, by Scalise, Olness and Tung [Scalise96]. 
In this approach the authors employ a strategy developed in the case of DIS in refs. 
[Aivazis94, Aivazis94a]. Here initial- and final-state mass singularities are resummed 
as in the fragmentation-function approach, and the result is then matched in the low- 
Pt region to the fixed-order NLO massive calculation. At large Pt this calculation 
does not include as yet, however, the full set of NLO corrections to the hard-process 
matrix elements. 

The preliminary numerical results of this study [Scalise96] are consistent with 
those of the approach by Cacciari and Greco; in particular, they support the conclu- 
sion that in the pr range explored by the Tevatron experiments the resummed results 
are consistent with the fixed-order ones, provided a scale n of the order of fio/2 or 
smaller is selected. 

The fragmentation function formalism discussed so far can also be applied to high 
transverse momentum charm- meson production [Cacciari96b]. 

4- 1-3 Comparison with experimental results 

The comparison of the final experimental data from UAl [Albajar91] with theory is 



shown in fig. ^ The data are (approximately) a factor of 2 higher than the central 
prediction. Good agreement can be obtained by selecting the largest theoretical 
prediction. 

Initially, the comparison with the Tevatron data led to pessimistic conclusions 
on the capability of NLO QCD to describe the data. Limited statistics in the chan- 
nels with fully reconstructed 6-meson decays [Abe92a], and the insufficient theoret- 
ical understanding of the J ftp production mechanisms, which was required at that 
time for the extraction of the b cross section from the measurement of inclusive 
J/ip^s [Abe92b, Abe93b], led to evidence that the production rate of b quarks exceeded 
by several times the one predicted by NLO QCD. The disagreement with theory was 
less pronounced in the case of higher-statistics measurements using semileptonic b 
decays [Abe93, Abe93a], which however probe production at pj, values higher than 
those based on exclusive and J/ip decays. 

The analysis of the large data samples accumulated during the Tevatron data 
taking of 1992-93 has recently been completed, and the first results from the data 
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Figure 19: Differential b-quark pj. distribution: comparison of the 
fixed- order NLO calculation (NDE stands for [Nason88]) with the frag- 
mentation-function approach [Cacciari94]; /^o is the scale at which the 
boundary conditions for the fragmentation functions are set. 



of the 1994-96 run are being released. This wealth of data and the improved de- 
tection capabilities offer today a rather precise picture of inclusive b production at 
1800 GeV. The CDF experiment primarily used its silicon vertex detector to improve 
the background rejection in the reconstruction of fully exclusive 6-hadron decays and 
to separate the J/t/^'s coming from b decays from those promptly produced, thus 
eliminating an important source of systematic error. The measurement of fully re- 
constructed 6-hadron decays also allows a precise determination of the differential 
Pt spectrum of b hadrons, free of the uncertainties related to the modelling of frag- 
mentation and decay. The experiment DO exploited its low muon backgrounds to 
extend the measurement of the b spectrum to high values of Pt [Abachi95a] , as well 
as using dimuon pairs [Abachi96, Abachi96a] and forward muons for a measurement 
of small-angle b production [Abachi96b] . 



The results obtained by the two collaborations are shown in figs. ^ and ^, 
compared with the theoretical predictions. 
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Figure 21: CDF data on the integrated b-quark pj- distribution, com- 
pared to the results of NLO QCD. 



For an easier evaluation of the results, we present them on a linear scale in the 



form Data/Theory in fig. ^, where we include the UAl data as well. We divided 
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Figure 22: DO data on the integrated b-quark pj. distribution, compared 
to the results of NLO QCD. 



the data by our central theoretical prediction. The dot-dashed straight lines are 
constant fits to the ratios, weighed by the inverse of the experimental uncertainties. 
For comparison, we also include the upper and lower theoretical curves divided by 
the central one (solid lines). 

The first important thing to notice is that, independently of the chosen input pa- 
rameters, the shape of the theoretical curves agrees very well with the data. Secondly, 
it must be pointed out that the results at yS = 630 GeV are consistent with those 
at vS^ = 1800 GeV. The average ratio Data/Theory measured by UAl differs by 
less than 10% from the ratio measured by DO, independently of the input parameters 
chosen. The difference is between 40 and 50% if one uses the CDF data. If we naively 
average the Tevatron data, we conclude that the relative discrepancy between theory 
and data when changing the value of v^ from 630 to 1800 GeV is about 30%, a 
number of the same order as the estimate of small-x effects. 

Independently of the beam energy, the data are higher by a factor of about 2 than 
the default prediction based on /i = /iq. They are, however, well reproduced by the 
upper theoretical curve. Therefore, while the overall uncertainty of the theoretical 
prediction due to the scale choice is large, there is currently no inconsistency between 
theory and data for the inclusive Pt distribution of b quarks at the Tevatron. The 
30% discrepancy between the results of CDF and DO is comparable to the discrepancy 
between the extrapolation of the UAl data to CDF, while UAl and DO data agree 
at the level of 10%. 



44 Heavy- Quark Production 



UAl < o > = 1.91 




10 20 30 

pr (GeV) 



60 



Figure 23: Linear comparison between experimental data and theory 
for the integrated b- quark pj. distribution. See the text for details. 



Recently, new data on b production at 630 GeV have been presented by CDF 
[Abe96a]. The measurement was performed during a special Tevatron run at Vo = 
630 GeV, using large-impact-parameter muons from 6-meson decays. The measure- 
ment was compared with a similar one carried out at y/S = 1800 GeV, to reduce 
common systematic errors. The preliminary result of this study is: 



abjPT > 9-5 GeV, |y| < 1, y^ = 630 GeV) 
ab{pT > 9.5 GeV, |y| < 1, v^ = 1800 GeV) 



0.193±0.025 Stat ± 0.023 syst . (4.5) 



To compare this result with a theoretical expectation, we have to choose the scales 
to be used at the two different energies. While there is no theorem on how to do 
this, it is however reasonable to choose the scales to be the same. This can be 
partly justified by the similarity of the scale dependence at the two energies, shown 
in fig. |18[ Assuming /i(630) = /i(1800) = /iq and using MRSA' parton densities, we 
find a theoretical ratio of 0.189 ± 0.012. The theoretical error is obtained by varying 
fj, in the range /io/2 < /i < 2/io and nib in the range 4.5 GeV < m^ < 5 GeV. 
The ratio (Data/Theory) (630)/(Data/Theory) (1800) measured by CDF is therefore 
equal to 1.0 ± 0.2. This is a very important measurement, because a large fraction 
of the systematic uncertainties present in the individual measurements cancel out in 
the ratio. This result indicates that possible sources of corrections to the theoretical 
calculations (e.g. higher-order corrections, small-x effects) should have the same 
impact at 630 and at 1800 GeV. One would reach the same conclusion by comparing 
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the results of UAl and DO, although the comparison of the systematic errors is in 
this case more complex. 

A recent measurement by DO [Abachi96b] , using muonic semileptonic decays of h 
quarks, explored the production of h quarks in the forward region (2.4 <\y^\ < 3.2). 
The results, shown in fig. |2^, indicate a rate in excess by a factor of approximately 4 
over the central theoretical prediction. The large samples collected at 1800 GeV have 
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Figure 24: Comparison between DO data and theory for the 

pp^{b^fi) + X cross section at large rapidity. 



also allowed a measurement of the differential Pt spectrum of fully reconstructed B 
mesons. CDF performed this measurement by detecting several hundreds of events 
in the final states B^J/^K^*^ [Abe92a, Abe95a, Abe96f]. The full reconstruction 
of the final state allows a precise measurement of the hadron momentum, therefore 
providing a direct measurement of the differential pj. spectrum. The systematic errors 
due to the modelling of the 6-hadron decays are significantly reduced relative to the 
measurements in the inclusive channels. The comparison with the NLO QCD calcu- 
lations, however, still requires a description of the 6-quark to 6-hadron fragmentation 
function. The NLO calculation of the 6-quark pj. distribution contains the effects of 
final-state gluon radiation, although limited to the emission of one gluon. As was 
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done in the fixed-target section, we then fold the NLO quark pj. spectrum with a phe- 
nomenological parametrization of the non-perturbative quark-to-hadron transition. 
We use the so-called Peterson fragmentation function [Peterson83] , already described 
in section 0. For the parameter e that defines the shape of the distribution we choose 
the value e^ = 0.006 ± 0.002 [Chrin87]. We furthermore assume that the fraction of b 
quarks that fragments into the hadron species identified in the experimental analyses 
(Bud) is 75%, a number supported by recent direct measurements [Abe96c]. 



Figure 25 shows the latest CDF data [Abe96f], compared to the NLO QCD predic- 
tions. The ratio Data/Theory is consistent with that found in the inclusive measure- 
ments, with the possible exception of the point at the lowest pr, which is higher than 
expected. It is premature to conclude whether this is just a statistical fluctuation, or 
whether this is a sign of small-x effects. 

As in the case of charm production at fixed target, we should comment here that 
the use of the Peterson fragmentation function might not be justified in the context 
of hadropro duct ion of heavy quarks. As a simple observation, we point out here that 
while the measurement of heavy-quark spectra in e~^e~ data is mostly sensitive to the 
first moment of the fragmentation function, corresponding to the average jet energy 
carried by the heavy hadron, the pr distributions in hadronic collisions are sensitive 
to higher moments of the non-perturbative fragmentation function. In fact, assuming 
for simplicity a perturbative pr spectrum of the form: 

apr Pt 

it is easy to prove that the resulting hadron spectrum, after convolution with a given 
fragmentation function f{z), will be given by: 

'^'^ 47 X [ dzz-'fiz). (4.7) 



n 



dpT Pt 

In the case of the Tevatron, n turns out to be approximately 5. It is not unlikely, 
therefore, that alternative models for the non-perturbative fragmentation of heavy 
quarks, which could equally well fit the e+e" data, could give rise to significant 
differences when applied to production in hadronic collisions. We also remark that 
the gluon component of the fragmentation function is not important in e~^e~ physics, 
while it may be crucial in hadropro duct ion. 

To display the size of the hadronization corrections induced by the Peterson frag- 
mentation, we also include in fig. ^the 6-quark p^ distribution. Notice that, neglect- 
ing the lowest-p-r point, both shape and normalization of the data are well described 
by the central theoretical prediction without fragmentation. 

We summarize the main conclusions of the studies presented in this section: 

1. There is good agreement between the shape of the 6-quark pj. distribution pre- 
dicted by NLO QCD and that observed in the data for central rapidities. 
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Figure 25: Comparison between CDF data and theory for the differ- 
ential B-meson Pt distribution. 



2. Although the data are higher by a factor of approximately 2 with respect to the 
theoretical prediction with the default choice of parameters, extreme (although 
acceptable) choices of A^^ and of renormalization and factorization scales bring 
the theory in perfect agreement with the data of UAl and DO, and within 30% 
of the CDF measurements. 

3. The choice of low values of /ipt and fip is favoured by studies of higher-order 
logarithmic corrections. 

4. There is a difference at the level of 30% between the measurements of CDF and 
DO at 1800 GeV, and between CDF and UAl at 630 GeV. 

5. The CDF measurements at 630 and 1800 GeV indicate that theory correctly 
predicts the scaling of the differential p^ distribution between 630 and 1800 GeV, 
a fact that had often been questioned in the past and now finds strong support. 

6. Forward production of b quarks indicates a larger discrepancy between theory 
and data. 

7. More theoretical studies should be devoted to the understanding of the non- 
perturbative fragmentation function for heavy quarks. As already stressed in 
the section on fixed target, there is no firm evidence, either theoretically or 
experimentally, that the standard Peterson parametrization is well suited for 
the description of the hadroproduction data. 
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4.2 66 correlations 

Studies of tlie one-particle inclusive distributions are a benchmark test of QCD. 
However, these tests do not probe all of the important features of the production 
mechanism. In several cases the correlations between the quark and the antiquark 
can provide additional important information. We already discussed how, in the case 
of fixed-target charm production, the shape of the azimuthal correlations between the 
c and c quarks is sensitive to the amount of intrinsic transverse momentum carried 
by the partons in the hadron. Given the large mass of the h quark and the large 
Pt values at which 6's are probed in hadronic collisions, these effects should be neg- 
ligible. However, one could observe similar behaviours induced by the build-up of 
large transverse momentum due to the perturbative evolution of the initial state (i.e. 
multi-gluon emission). This phenomenon would be even more pronounced at small 
Pt and in the presence of important small-x effects. 

The knowledge of heavy-quark correlations is also an important element in the 
study of detector acceptances and detection systematics. For example, in studies 
of h tagging in samples where both 6's have been tagged, it is important to know 
what the correlations between the h and h momenta are, since the tagging efficiency 
depends on the h momentum. These studies have implications on the determination 
of tagging efficiencies for a multitude of important measurements, from the study of 
top production [Abe94] to the study of R^ [Nason96] . 

The study of hb correlations in hadronic collisions was pioneered by UAl [Geiser92, 
Albajar94, Albajar96]. For these measurements they used their sample of high-mass 
and non-isolated dimuon events. The shapes of the A0 and Ai? distribution for the 
data are both in good agreement with the theoretical expectations [Albajar94]. The 
rate is approximately 30% higher than the central value of the theoretical predictions, 
a result consistent with the UAl determination of the inclusive cross section. They 
also separated the sample of dimuon events into two classes, defined by two-body and 
three-body topologies. Contributions to the three-body topology come from events 
with hard-gluon bremstrahlung, and are therefore proportional to a|. Two-body final 
states, with the h and the h almost back-to-back, are dominated by Born-like processes, 
and their rate is of 0{a1). A comparison between the two samples, carried out as a 
function of the Pt of the h quark, allowed UAl to study the dynamical evolution of hb 
correlations and to extract the value of as{Mz) = 0.113]^g[jgg ^0009 i, [Albajar96]. 
The same study led to a 3o" evidence for the running of as- 

CDF presented their first study of hb correlations using a sample of e/i pairs 
[Abe94a]. Once again, good agreement was found between the shape of the data 
and the QCD expectations. A second measurement has been recently reported by 
CDF, using muons plus tagged b jets [Abe96b]. The shape of the A0 distribution is 
in reasonable agreement with QCD, while the distributions of the jet E^^ and of the 
muon Pt are not (see fig. pG]) . 

Contrary to the case of cc correlations measured in fixed-target experiments, the 
measurement at hadron colliders is sensitive to the modelling of the heavy-quark 



Hadron colliders 49 



:dCT/dET 

2 - 



M^ \ 



;8 o o o 



-X) er- 



Jet Et (GeV) 



O 20 



40 50 60 70 

Pt (GeV/c) 



(a) 



10 20 30 40 50 60 70 80 90 100 




90 100 

















1 


= da /66(p 


* ' 


f + 


* 


« 


(c) 




-- O I) 

T 


■ o 


o 





• 


1 


l1 










o 


o 


- 1 1 1 1 1 1 


1 1 1 1 


1 1 1 1 


, 


, 


,1 



0.5 



Radians 



Figure 26: CDF results on the bb correlations using ^+b-tagged jet 
final states [Abe96b], compared to NLO QCD. Solid points correspond 
to the default theoretical prediction, with scales fi = fio, empty circles 
correspond to the difference between the choice fi = fJ,o/2 and /x = /xq. 
Top figure: tagged-jet Et distribution. Central: muon px spectrum. 
Bottom: azimuthal correlations. 



fragmentation, because of possible trigger biases. A harder (softer) fragmentation 
function would enhance (decrease) the efficiency for the detection of the softest of 
the two b quarks. These effects could have an impact on the distributions reported 
in this study by CDF. This collaboration explored the effects of changes in the e?, 
parameter within the Peterson fragmentation model, finding them negligible. As 
we argued earlier, it cannot be excluded that a systematic study of other possible 
parametrizations for the fragmentation modelling could lead to significant effects. 
Also the possible effects of the kr kick have been studied by CDF [Abe96b] , with the 
conclusion that not even an average kj, as large as 4 GeV could improve the agreement 
between theory and data for the measurements considered. 

Both CDF and DO recently presented studies based on samples of high-mass 
dimuons [Abe96d, Abachi96a]. The shapes for both A0 and p^ for a given p^ are well 
reproduced by theory, within the large uncertainties, while there is a normalization 
discrepancy relative to the central predictions (fig. 
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Figure 27: Results on hh correlations using diniuon final states, com- 
pared to NLO QCD. Left: azimuthal- correlation results from DO. Right: 
Pj. distribution of the slowest muon from CDF. 



In conclusion we can say that, overall, there are good indications that NLO QCD 
correctly describes the correlations between b pairs produced in hadronic collisions. 
Most of the large theoretical systematics present in the case of fixed-target charm 
production are not relevant here. The discrepancies observed in the CDF measure- 
ment using the muon+tagged-jet final states, however, indicate that there still is a 
possible interplay between the theoretical and experimental systematics, which is not 
entirely understood as yet. 



4.3 Heavy-quark jets in perturbative QCD 

4.3.1 Preliminaries 

An interesting way of understanding the production mechanism of heavy flavours is 
to consider the cross section of jets that contain a heavy quark (in short: heavy-quark 
jets) [Frixione96] . The main difference between the study of a heavy quark and of a 
heavy-quark jet is that in the former we are interested in the momentum of the quark 
itself, regardless of the properties of the event in which the quark in embedded, while 
in the latter we are interested in the properties of a jet containing one or more heavy 
quarks, regardless of the momentum fraction of the jet carried by the quark. A priori 
it is expected that variables such as the Ej, distribution of a heavy-quark jet should be 
described by a finite-order QCD calculation more precisely than the pj. distribution 
of open quarks. This is because the large logarithms log{prp/mh) arising from the 
emission of collinear gluons from the heavy quark cancel in this case. Furthermore, 
the experimental measurement of the Ej. distribution of heavy-quark jets does not 
depend on the knowledge of the heavy-quark fragmentation functions, contrary to the 
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case of the pr distribution of open heavy quarks. Experimental systematics, such as 
the knowledge of decay branching ratios for heavy hadrons or of their decay spectra, 
are also largely reduced. 

The calculation of the heavy-quark jet rate is very similar to the one of the generic 
jet cross section. Two important differences have nevertheless to be stressed: by its 
very definition, a heavy-quark jet is not flavour-blind; one has to look for those jets 
containing a heavy flavour. Furthermore, the mass of the heavy flavour is acting 
as a cutoff against final-state collinear divergences. This in turn implies that the 
structure of the singularities of the heavy-quark jet cross section is identical to that 
of the open-heavy-quark cross section. The heavy-quark jet cross section at NLO can 
therefore be written in the following way, which implicitly defines dA/dEj, (this term 
will be denoted in the following as the "jet-like" component of the cross section) 

da rfCT(°P^") dA ,, , 

(4.8) 



dEj< dErp dEj. 

where da^°^^'^' / dE^ is the open-heavy-quark cross section. Observe that dA/dE^ 
has no contributions at the Born level. Gluon-radiation effects start contributing to 
it at order q;|. Normally the heavy-quark jet is defined by a cone algorithm. In 
our phenomenological applications, we will use the Snowmass convention [HuthQl], 
whereby particles are clustered in cones of radius R in the pseudorapidity-azimuthal 
angle plane. One then requires that the h quark be inside the jet. At the NLO, the 
heavy-quark jet can be the heavy quark, or it can contain the heavy quark and the 
light parton, or the heavy quark and the heavy antiquark. In other words, there is 
the possibility to cluster more than one parton to get the heavy-quark jet, which 
will be eventually observed. It turns out that, if the jet definition is appropriately 
infrared-safe, all the subtractions needed to get an infrared-safe result are contained 
in the term da'^"'^^'^^ /dEj.. This is shown in details in ref. [Frixione96] . This result is 
however quite intuitive. In fact, singularities arise because of radiation collinear to 
the incoming partons, which cannot enter the jet cone, or to soft radiation, which 
does not sensibly alter the jet Et- 

4-3.2 The structure of heavy- quark jets at the Tevatron 

We present in this section some results of interest for the Tevatron collider [Koehn95] . 
We consider jets produced within |?7| < 1, in order to simulate a realistic geometri- 
cal acceptance of the Tevatron detectors. We will use the parton distribution set 
MRSA' [Martin95a]. Our default values of the parameters entering the calculations 
are mc = 1.5 GeV, m^ = 4.75 GeV and /x^ = /ij? = /iq = JEt + "^q- I^ practice, since 
we will consider jets of energy much larger than the heavy-quark mass, the depen- 
dence of our results upon the mass value is almost negligible. We will therefore study 
only the uncertainties associated with the choice of renormalization and factorization 
scales. 
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Figure 28: b jet inclusive Erp distribution in pp collisions at 1.8 TeV, 
for \ri\ < 1, R = 0.7 and /ip = fi^ = fiQ (solid line). For comparison, we 
also show the open-quark inclusive Ej. distribution (dashed line). The 
component of the jet-like contribution due to jets containing both b and 
b is represented by the dotted line. 



Figure ^ shows the prediction for the Et distribution of b jets at the Tevatron 
for R = 0.7. For the purpose of illustration, the open-quark component is presented 
separately. It is apparent that the jet-like component becomes dominant as soon as 
Erp becomes larger than 50 GeV. It can be shown [Frixione96] that this value depends 
significantly on the cone size, being equal to 25 and 100 GeV for R = 1 and 0.4 
respectively. Also, for a fixed R, the value of E^ at which the jet-like component 
becomes dominant is smaller for c jets than for b jets. This is in agreement with naive 
expectations. The relative probability of finding a heavy-quark pair inside a high-Erp 
gluon scales in fact like \og{Erp/mQ) [Mueller85, Mueller86, Mangano92, Seymour95]. 

We also show the part of the jet-like component due to jets including the bb pair 
(we will call these bb jets). The figure suggests that for this Ej. range and with 
R = 0.7 this is the dominant part of the jet-like component. This is consistent with 
the expectation that, for large enough Ej., and provided that the majority of the final- 
state generic jets are composed of primary gluons, heavy-quark jets are dominated 
by the process of gluon splitting, with the jet formed by the heavy-quark pair. As we 
will show later, the situation changes at higher values of Ej, where heavy quarks are 
mostly produced via the s-channel annihilation of light quarks. 



The left-hand side of fig. |29| presents the theoretical prediction for the absolute b- 
jet rate at Ej. = 50 GeV versus the cone size, for different choices of the factorization 
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Figure 29: Left: b jet inclusive Ej. rate, as a function of the cone size 
R, at Erp = 50 GeV and for various scale choices (fi^ = fip = jjl). Right: 
Scale dependence of the b-jet Ej. distribution (R = 0.4, solid lines) and 
of the open- quark inclusive E^ distribution (dashed lines). 



and renormalization scales. The cross section at i? = is well defined, and it is 
equal to the open-bottom cross section. This should be contrasted with the case of 
generic jets, in which the cross section at i? = is not well defined, being negative 
at any fixed order in perturbation theory [Aversa89, AversaQOa, Ellis89S, EUisQOS]. 
Completely analogous results hold in the case of c-jet production. The right-hand 
side of fig. ^ shows the scale dependence of the 6-jet cross section (i? = 0.4) as a 
function of Et, for values up to 450 GeV. 

The strong scale dependence exhibited by the absolute rate at low and moderate 
i^T values is of the same size as the one present in the inclusive p^ distribution of 
open bottom quarks. This scale dependence is usually attributed to the importance 
of the gluon splitting contribution. One therefore expects that, in a regime where 
the gluon-splitting contribution is suppressed by the dynamics, the scale dependence 
should be milder. We will show later that this suppression does indeed take place 
for high transverse energies. This explains why, in the high-£^T region, the scale 
dependence is indeed reduced to the value of 20% when the scale is varied in the 
range /io/2 < yU < 2/io, a result consistent with the hmited scale dependence of the 
NLO inclusive-jet cross sections [Aversa89, AversaQOa, Ellis89S, Ellis90S]. In spite 
of the strong scale sensitivity at the smaller values of -Etj it is reasonable to expect 
that the ratio of the h- and c-jet rates be a stable quantity under scale variations. It 
can be shown [Frixione96] that this is indeed the case: by taking /io/2 < ^ < 2/io the 
ratio changes at most by 10% with respect to the default prediction; furthermore, the 
ratio of 6 to c jets is more stable than the one of 6 to c quarks. 

Of direct interest for studies of heavy-flavour tagging and for searches of possible 
new physics is the fraction of heavy-quark jets relative to generic jets. This is also 
in principle the most straightforward measurement from the experimental point of 
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Figure 30: Ratio of the heavy-quark jet (for bottom, left, and charm, 
right) to inclusive-jet Ej. distributions, for different choices of renor- 
malization and factorization scales (fin = fip = fi), for R = 0.4 (top) 
and R = 0.7 (bottom). The data points for R = 0.4 represent prelim- 
inary results from the CDF experiment, for which only the statistical 
uncertainty is shown. 



view. We present in fig. ^ the ratio of the heavy-quark jet to the inclusive-jet Ej, 
distributions [Aversa89, AversaQOa, Elhs89S, ElhsOOS] for bottom (left) and charm 
(right). The inclusive-jet E^ cross section used here was calculated with the JETRAD 
program [Giele94], using the same choices of parton densities and {fifi,fip) that were 
adopted for the 6-jet and c-jet calculations. Contrary to the figures presented so far, 
which showed results for the heavy quark only (i.e. no antiquark contribution), we 
adopt for this figure the prescription used in the definition of the data presented by 
CDF [Koehn95]. The Q-jets are defined there as jets containing either a Q ot a Q 
quark, jets containing both being counted only once. We will call these Q{Q) jets. 
This distribution can be obtained by subtracting the contribution of QQ jets from 
twice the total Q-jet rate. 

It is interesting to notice that, as far as b production is concerned, there is a 
good agreement between the CDF data and the theoretical prediction obtained with 
fin = ftp = /io/2; this choice of scale is also supported by the inclusive-jet i?T-spectrum 
data [Abe92, Abe96, Blazey96]. This is particularly significant since the choice of scale 
for the heavy-quark jet cross section is not independent from the scale chosen to pre- 
dict the open-heavy-quark one (see eq. ( ^^) ). On the other hand, to get a satisfactory 
description of the data for the open-bottom pj. spectrum, smaller scale values (of the 
order of /io/4) or larger Aqcd values (compatible with LEP measurements) have to 
be chosen [Frixione94, Frixione96a]. Should this situation persist when additional 
data on b jets will become available, it would indicate an inconsistency in describing 
two phenomena due to the same underlying physics. The poor understanding of the 
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Figure 31: Ratio of the h{b) jet to inclusive-jet Ej- distributions for 
flu = fip = fj,/2 and with R = 0.4 (solid) or R = 0.7 (dashes). The data 
points are preliminary CDF data [Koehn95j, obtained with R = 0.4. 



fragmentation mechanism is very likely a source of this inconsistency. As for the 
large disagreement with the charm data, we have no significant comment to make. 
Hopefully, additional data will soon be available, as well as estimates of the experi- 
mental systematics. Notice that the largest contribution expected from higher-order 
perturbative corrections is given by the production of cc pairs from soft gluons emit- 
ted during the gluon-shower evolution. However, these effects have been estimated 
in refs. [Mueller85, Mueller86, Mangano92, Seymour95], and have been shown to be 
negligible at the energies of interest for the current measurements. 

To conclude this section, we discuss the behaviour of the 6-jet production cross 
section at high Et- The interest of this item stems from the discrepancy reported 
by CDF [Abe92, Abe96] in the tail of the jet distribution. If this discrepancy could 
not be accommodated by new theoretical developments in QCD [Catani96a] or in the 
fitting of parton densities [Huston95, Glover96], a study of the flavour composition of 
these high-energy jets could help in understanding the nature of the phenomenon. 



Figure ^ shows the 6(6)-jet fraction for Ej, values up to 450 GeV, for two different 
values of the cone size {R = 0.4 and R = 0.7) and at fia = fip = /io/2. Notice 
that while the fraction remains constant through most of the Et range, a rise is 
observed above 300 GeV. To better understand the origin of this rise, we present 



in fig. p^ a the separate contribution to the 6-jet cross section of the three possible 
initial states, gg, qq and qg. Notice that the qq contribution becomes dominant 
for Ej. > 250 GeV. Figures [5^-d show, for each individual channel, the separate 
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Figure 32: Initial-state composition of tfie b-jet production processes, 
calculated for ^u^ = fip = fio/2 and R = 0.4 (upper left). Different com- 
ponents of the production processes: gg -^ b-jet (upper right), qq ^ b 
jet (lower left) and qg ^ b jet (lower right). 



contribution of the open-quark and 66-jet components. For Et large enough, the 
dominant component of the gg and qg channels is given by the 66-jet contribution, 
because of the gluon-splitting dominance. In the case of the qq channel, on the 
contrary, the bb-jet term is always suppressed, and most of the b jets are composed of 
a single b quark, often accompanied by a nearby gluon. We conclude that at high E^ 
the dominant mechanism for the production of a 6 jet is the s-channel annihilation of 
light quarks. Since mass effects are negligible at high Et, 1/5 of the jets produced in s- 
channel annihilation are 6 jets. A simple LO calculation shows that the fraction of the 
two-jet rate due to s-channel light-quark annihilation is about 20% at Et = 450 GeV, 
giving an overall 6-jet over inclusive-jet fraction of approximately 4%. This explains 
the rise of the b fraction at high Et, and provides a nice consistency check of our 
results. Notice also that while the probability that a gluon jet will split into a bb pair 
grows at large Et faster than the 0{a^) result [Mueller85, Mueller86, Mangano92, 
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Seymour 9 S]!^! the fraction of primary gluons in the final state is so small that the 
overall effect on our predictions is negligible. 

4.4 Associated production of heavy quarks with W or 7. 

In addition to the standard measurements of inclusive heavy-quark production, the 
high energies and luminosities available at the Tevatron collider allow the detection 
of the associated production of heavy quarks and electroweak vector bosons (7, W 
and Z). The study of these processes goes beyond the goal of testing perturbative 
QCD, and has a wide range of applications, which will be briefly described in this 
section. 

4-4- 1 Photon plus heavy quarks 

The associate production of a photon recoiling against a jet is dominated by the 
Compton scattering of a gluon and a quark in the proton sea. This process, in 
the most common experimental configurations, dominates over the quark-ant iquark 
annihilation channel. In the case of a heavy-flavoured jet we would then be sensitive to 
the heavy-flavour component of the proton sea. The measurement of the Pt spectrum 
of photons produced in association with charm jets, and the knowledge of the gluon 
parton density, could therefore provide in principle a direct measurement of the charm 
density inside the proton [Fletcher89] . The similar process gb^'-fb would allow the 
study of the b density. In addition, these processes provide a useful control sample 
for the study of heavy-quark tagging. 

Figure |33| contains the distributions of jets of various flavours produced with 
photons, as estimated with a LO QCD calculation and CTEQIM parton densi- 
ties [Botts93]. Because of the difference in charge and partonic density, associate 
production with a bottom quark is suppressed by a factor of 8-10 relative to charm. 
Notice also the curious fact that, because of the suppression of the light-quark anni- 
hilation channel, it is more likely for a jet produced in association with a photon to 
be a charm jet rather than a gluon jet, at least for transverse momenta up to 30 GeV. 

There is a theoretical interest in these processes because the density of a heavy 
quark inside the proton is in principle calculable perturbatively. Neglecting higher- 
order logarithmic corrections, which can be resummed using the Altarelli-Parisi evo- 
lution, the inclusive process pp—>-'-fQ + X can be calculated by evaluating the par- 
tonic process gg^QQ'~f and integrating over the phase space of the Q [Mangano94, 
Stratmann95]. This process is dominated by configurations where the quark being 
integrated over is produced at large rapidity and small Pt- No divergence will appear, 
because of the heavy-quark mass. A consistent definition of the heavy-quark density, 
including threshold effects [Collins86] , should then reproduce the result of this calcu- 
lation, at least at moderate Pt- Comparison against the experimental data is however 
important, to verify that no additional non-perturbative effect is at work. 
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Figure 34: Photon p^ distribution in events with a central b jet in 
direct photon events, showing the results of the LO structure-function 
approach and of the 0{a'^) calculation. 
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We show in fig. ^ the photon Pt distribution in events with a central 6, calcu- 
lated [Mangano94] using the structure function approach {bg^b'j, sohd line, CTEQIM 
PDF's [Botts93]) and using the exact matrix elements for the gg + qq^'yQQ pro- 
cesses [Ellis88] (dashed line). The qq channel produces a heavy-quark pair via gluon 
splitting, and cannot be accounted for by the structure-function calculation. As the 
plot shows, there is perfect agreement between the two approaches, at least in the 
region where gluon splitting is small. This indicates that the effects of initial-state 
evolution for the b quark at these values of x and Q^ are not important. Notice that 
this conclusion is important for the consistency of the NLO evaluation of the inclusive 
b cross sections, which contains only the flavour excitation diagrams present in the 
0{a^) matrix elements. 

Similar results are obtained in the case of charm production, where there is, how- 
ever, some larger sensitivity to the choice of the charm mass. A more comprehensive 
study of the charm case can be found in ref . [Stratmann95] . 

A NLO QCD calculation of this process has also been recently performed [Bailey96]. 
This calculation uses massless charm quarks; it is therefore valid in the limit of large 
transverse momentum. The inclusion of photon isolation and of a c^7 fragmenta- 
tion function is then necessary to remove the divergences due to the coUinear emis- 
sion of a photon from the final-state charm quark. The first experimental results 
on the associated production of photons and charm have recently appeared from 
CDF [Abe96e]. The quoted result is the cross section for the process pp-^jD*^ + X, 
for Pt(7) > 16 GeV, Pt{D*+) > 6 GeV, \r]{-f)\ < 0.9 and \ri{D*+)\ < 1-2. The re- 
ported measurement is a = 0.38 ± 0.15 ± 0.11 nb. While the statistics of this study 
(based on approximately 19 pb~^, i.e. a 20% fraction of the full data set currently 
available) are not sufficient to produce a distribution in Pt, the agreement of theory 
and data is satisfactory. Estimates performed using the LO PYTHIA Monte Carlo 
program [Sjostrand87], using as renormalization scale /i = Pri'j), give rates in the 
range 0.18-0.22 nb, depending on the set of parton densities used. Given the large K- 
factor values found in ref. [Bailey96] between the LO and NLO predictions, the CDF 
results disfavour a significant presence of a non-perturbative component of charm 
quarks inside the proton in the regions of x explored by this measurement. 

4 -4 -2 W bosons plus heavy quarks 

The associated production of W^ bosons and heavy quarks is interesting for several 
reasons. Production of W plus charm is, in analogy to the 7 plus charm case, a 
potential probe of the partonic densities of the proton. In this case, the Cabibbo- 
allowed channel gs^Wc is the leading production mechanism, and provides a direct 
probe of the strange content of the proton at large Q^ [Baur93]. 

The Wc final state is also one of the dominant sources of events with a W plus 
one jet with a reconstructed secondary vertex. As such, these events enter in the 
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This happens because of pairs emitted at higher orders from the gluons of the shower. 
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analysis of the low jet multiplicity W^+jet events, which is performed as ancillary 
work to the isolation of the ti signal (see the extensive discussion in ref. [Abe94]). 
Although no cross sections for the W + c process have been explicitly quoted by the 
experiments, the agreement of the data with the LO theoretical predictions, modelled 
for the presence of several analysis cuts, is good [Abe94, Abe95, Abachi95]. The full 
set of NLO (0(a^)) corrections to this process has recently been evaluated [Giele96], 
resulting in an improved stability with respect to renormalization and factorization 
scale variations. 

Because of the small value of Vcb and of the charm density in the proton, which 
suppress the LO process gc^Wb, associate production of W and bottom quarks 
has its dominant contribution from the higher-order process qq—^Whh [Kunszt84, 
Mangano93a]. Interest in this process is once again related to its relevance for the 
isolation of the ti final state, which has always among its decay products a W boson 
and a hh pair. Direct evidence for the prompt production of Whh events comes 
from the study of tagged jets in the W + 2J sample, performed by the CDF and 
DO collaborations in the top analyses (for a detailed discussion, see ref. [Abe94]). 
Because of the small tagging efficiencies for charm quarks, and because of the limited 
statistics, there is currently no experimental result on the production of Wcc events. 

Additional interest in the Wbb process comes from the studies of direct back- 
grounds to the detection of an intermediate-mass Higgs. This could be observed with 
possible future luminosity upgrades of the Tevatron collider in the mass range 65 < 
Mfi < 110 GeV using the pp-^WH-^Wbb channel [Stange94, Stange94a, Amidei96]. 
The DO collaboration has recently presented the first study of a direct search for 
Higgs bosons in the Wbb channel [Abachi96d]. The findings are consistent with the 
QCD background estimates, providing preliminary evidence for the correctness of the 
QCD calculations of the Wbb rate, although no explicit cross section has as yet been 
quoted for this process. 

4.5 Production of top quarks 

4-5.1 Total ti production cross sections 

The top quark having been found [Abe94, Abe95, Abachi95], the comparison between 
its observed production properties and those expected from the Standard Model will 
be an important probe of the possible existence of new phenomena. One of the most 
important tests to be performed concerns the total production cross section. This is 
the most inclusive quantity available, and is a priori the least sensitive to a detailed 
understanding of the higher-order corrections that infiuence the evolution of the initial 
and final states. Within QCD, one expects the perturbative expansion in powers of 
as{mt) to be well behaved and to provide an accurate estimate of the total cross 
section already at low orders. In particular, the first estimates of the total production 
cross section using the full NLO matrix elements [Nason88, Beenakker89, Ellis91] gave 
an increase (relative to the Born result) of the order of 30% for masses above 100 GeV. 
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The residual perturbative QCD uncertainty, evaluated by varying the renormalization 
and factorization scales, was estimated to be no larger than 10%. The choice of 
parametrization for the input parton densities was also shown to give effects of this 
order of magnitude, by using the available sets. 

It was later pointed out in refs. [Laenen92, Laenen94] that logarithmic contribu- 
tions associated with the emission of soft gluons from the initial state could signifi- 
cantly enhance the NLO result. Similar conclusions were reached in refs. [Berger95, 
Berger96]. More recent studies [Catani96, Catani96a], which will be described in more 
detail in Section ^, have proved that the effect of soft-gluon resummation is actually 
very small, of the order of a 1% correction to the total cross section. 

Beside the soft-gluon emission effects, also Coulomb effects [Sommerfeld39] may 
enhance or deplete the cross section near threshold, as originally discussed in ref. 
[Fadin90]. These effects can be evaluated in the following way. We separate the 
partonic Born cross-section formulae into their colour singlet and octet components 



'(8) ^"(1) 
'(8) 

where a^'^^'^^^ can be found in ref. [Nason88], and 



cr(««") = afsf (4.10) 



^(99) (^-s _ «i /^P^ 



1, 1 + /3, o, 



(4.11) 



(1) "'^ ~ m2 384 
Here (3 = y/1 — p and p = Am^/s. The Coulomb-resummed cross section is given as 

^ (P) = ^(8) (P) 7 TT^^^TT T + (^{1) (P) -. ^-. TT^^^ • (4. 12) 

exp(7ras/(6p)) — 1 1 — exp(— 47ras/(3/5)) 

We do not include bound-state effects, which, as shown in ref. [Fadin90], are much 
smaller. As shown in ref. [Catani96], the relative correction due to the resummation 
of Coulomb effects not already included into the NLO results is an effect smaller than 
1% of the NLO cross section. Together with soft-gluon resummation effects, higher- 
order Coulomb corrections can therefore be entirely neglected, given the theoretical 
NLO uncertainties and the current experimental accuracy. 

Electroweak processes can give rise to single top production. This happens through 
the so-called Wg fusion process [Dawson85, Willenbrock86, Yuan90, Ellis92, Bordes95], 
whereby a virtual W, emitted from an initial-state light quark, interacts with an 
initial-state gluon and produces a tb pair, and through the decay of an off-shell W 
boson into a tb pair [Cortese91, Stelzer95]. While in principle these processes do not 
contribute to the observation of ti pairs, and therefore should not interfere with the 
measurement and study of the QCD production properties of top quarks, in practice 
they constitute an irreducible background to the detection of ti events. From the 
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Figure 35: Scale dependence of the top cross section at tfie Teva- 
tron, for mt = 175 GeV. Dotted lines are the /Xp dependence at fixed 
jiji, dashed lines are the fj,^ dependence at fixed fip, and solid lines are 
obtained by the simultaneous variation of fj,^ and /ip. 



point of view of QCD studies, it is therefore a welcome fact that their total contri- 
bution to the inclusive top-production cross section is only a fraction of the order of 
25%. Their study is however important for electroweak measurements, as it probes 
the Wtb vertex and the value of Vtb directly. A complete review of the phenomenolog- 
ical applications of single top production, in addition to a complete list of the relevant 
references, can be found in ref. [Amidei96]. 

Electroweak radiative corrections to top production in hadronic collisions have 
been considered in refs. [Kao93, Stange93, Beenakker94, Harlander95] . They range 
from —0.97% to —1.74% of the Born cross section, for a Higgs mass of 60 and 
1000 GeV respectively^- 

In this section we will confine our theoretical analysis to results obtained within 
the QCD 0{al) approximation. It should be pointed out that this study is performed 
within the strict domain of the Standard Model. Corrections to the top cross section 
much larger than the Standard Model QCD uncertainties can be obtained in the pres- 
ence of new phenomena. Virtual corrections due to loops of supersymmetric particles 
have been studied, for example, in refs. [Li95, Li96, Alam96, Kim96, Sullivan96]. For 
a partial list of examples of alternative production mechanisms for top quarks, arising 
in theories beyond the SM such as Supersymmetry or Technicolor, see for instance 
refs. [Eichten94, Hill94, Holdom95, Casalbuoni96, Kane96]. 

^"'^W. HoUik and D. Wackeroth, private communication. 
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Figure 36: Top cross section as a function of as{Mz)- The dashed 
line (obtained with the MRS A' set) does not include the variation of the 
parton densities due to the change in as. 



We now proceed to a review of the theoretical uncertainties on the total cross sec- 
tion for tt production, which arise at NLO. Uncertainties due to unknown higher-order 
effects are usually accounted for by varying the renormalization (/Xj^) and factorization 
(/ip) scales. In principle, independent variations of the two scales should be consid- 
ered. In fig. ^ we show the scale dependence of the top cross section. Notice that 
there is a compensation of the renormalization-scale dependence when the different 
subprocesses are added up. In fact, the renormalization-scale dependence in the gg 
and qq channels has a behaviour opposite to that in the qg channel. It is only the 
combined scale dependence that can be considered an estimate of the neglected sub- 
leading corrections. As a second point, we observe that the maximum of the cross 
section is reached around mj/2. Thus, the usual choice of the range mt/2 < fi < 2mt 
appears to be particularly justified in this case. As a third point, we observe that the 
scale dependence in the cross section goes in the same direction for the two scales, 
so that, for the purpose of estimating the associated uncertainty, it is sufficient to 
consider the simultaneous variation of /Ir and /ip. 

Aside from the scale uncertainties, which reflect the limitations of the perturbative 
QCD calculation, there are uncertainties associated with our imprecise knowledge 
of the physical parameters involved. In particular, the strong coupling constant is 
determined within a certain accuracy. 

In the case of top production at the Tevatron, there is fortunately a compensating 
mechanism that reduces the dependence of the cross section upon as- In fact, the 
top cross section is dominated by the qq annihilation process. Quark densities in 
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4.13 


4.27 
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4.52 



Table 3: Total cross sections (in pb) for mt = 175 GeV at NLO. The 
set CTEQ' is taken from ref. [Huston95]. 



the proton are directly measured in DIS at a scale around 10 GeV^, and the QCD 
evolution to scales of the order of the top mass makes them softer. Therefore, for larger 
coupling, the partonic cross section increases, but the quark luminosity decreases. 

In order to perform a fair estimate of the uncertainty due to Aqcd, we need sets 
of parton densities fitted with different values of the strong coupling. The sets of 
ref. [Martin95a] meet our purpose^. In fig. ^ we show the cross section as a function 
of the strong coupling. 

For comparison we also show the ag dependence if the parton densities are kept 
fixed to the MRSA' set. We see the remarkable reduction in the as dependence, due 
to the compensation of the rise of the partonic cross section and the decrease of the 
quark parton densities. 

Our results for the top cross section for mt = 175 GeV are collected in table |. We 
also show results obtained with the recent parametrization derived in ref. [Huston95] 
by including the CDF jet data. We stress that the numbers in the table are obtained 
with a standard one-loop calculation. No resummation effects are included. 

Our range for the top cross section at mt = 175 GeV {mt = 170 GeV) is given by 



4.75 



+0.73 



7+0.86 



-0.62 P^ (^•^'''-0.73 

results: 



pb). This should be compared with the current experimental 



(Tttijnt 
(Ttt{mt 



175 GeV) 
170 GeV) 



7.5 ± 1.6 pb (CDF [Caner96]) 
5.2 ± 1.8 pb (DO [Narain96]) . 



(4.13) 
(4.14) 



^^The values of A4 that accompany the FORTRAN program for the structure function sets of 
ref. [Martin95a] are not consistent with the values of as quoted there, the differences being of the 
order of 1%. In the present work, we extract the values of A5 from their quoted values of as using 
the standard two- loop formula [Barnett96] . 
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Figure 37: Top cross-section at the Tevatron at \/5 = 1.8 TeV. The 
solid line is obtained with MRSA' parton density, and the dashed lines 
correspond to the upper and lower values obtained in table |j. The ex- 
perimental data are taken from refs. [Caner96, Narain96j. 



As a central value for our determination we have chosen the MRSA' [Martin95] result 
with yUR = /ip = rrit, in association with a value of A5 = 0.152 GeV (which corre- 
sponds to as{Mz) = 0.1113, according to the standard two-loop formula [Barnett96]). 
For the MRS sets with variable A, we have used A5 = 0.0994, 0.140, 0.190, 0.253, 
0.328 GeV (which correspond to as{Mz) = 0.105, 0.110, 0.115, 0.120, 0.125). 

The cross section bands are also shown in fig. 37. The agreement of theory and 
data is good, but it is clear that higher statistics should be collected before a significant 
test is achieved. 



get a{tt) 



For reference, we also quote the cross section for top production at the LHC. We 
i.77|'|g-^2 nb, for rrit = 175 GeV and v^ = 14 TeV. The error is obtained 
with the same scale and as variations that we used for the Tevatron cross section. 
Notice that the relative error is much larger than at the Tevatron, because of the 
greater importance of the gg initial-state contribution at the LHC energy. 



4-5.2 Top kinematical distributions 

Because of the high precision of the theoretical prediction, the measurement of the 
top cross section could become a sensitive probe of new physics. The only way 
to unambiguously disentangle the standard sources of top production from possible 
new physics will nevertheless rest on the detection of significant discrepancies in 
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Figure 38: Invariant-mass distribution of the tt pair. The HERWIG 
prediction has been rescaled by a constant factor K = 1.34. 



the kinematical features of top production from what is expected from QCD. For 
example, the detection of mass peaks in the invariant mass distribution of top pairs 
would indicate the possible presence of s-channel resonances strongly coupled to the 
top [Eichten94, Hill94]. 

Kinematical distributions for top-quark pairs produced at the Tevatron have been 
calculated at NLO in ref. [Frixione95] . Since top quarks manifest themselves ex- 
erimentally in a rather indirect way, their identification relies on complex series 
of experimental cuts [Abe94, Abe95, Abachi95, Campagnari96, Sinervo96]. The 
study of their kinematical properties therefore strongly relies on the modelling of 
the top-quark production mechanism itself. In particular, experimental analyses 
are performed using parton shower Monte Carlo (MC) programs, such as HER- 
WIG [Marchesini88, Marchesini92], ISAJET [Paige86] or PYTHIA [Sjostrand87]. It is 
therefore important to compare the distributions predicted by these MC calculations 
with what is expected from NLO QCD. The study performed in ref. [Frixione95] indi- 
cates an excellent agreement between these results, at least for distributions that are 
non-trivial at LO in perturbation theory, such as the top inclusive Pt distribution or 
the invariant-mass distribution of the top-quark pair (fig. ^). The theoretical uncer- 
tainty from higher-order corrections and from hadronization effects is very small for 
these variables, making them a rather solid term of comparison in the search for new 
physics. Distributions that are non-trivial only at (9(ag), such as the p^ of the top 
pair, show instead significant higher-order corrections. This is visible, for the example 
of the top-pair pj. distribution, in fig. ^, where a comparison is made between the 
NLO QCD predictions and those of the HERWIG MC. The multiple gluon emission 
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Figure 39: Transverse-momentum distribution of the tt pair. The 
HERWIG prediction has been rescaled by a constant factor K = 1.34. 



from the initial state of the hard process, included in the MC calculation, significantly 
increases the average transverse momentum of the pair, as can also be seen by the 
numbers quoted in table ^. Notice the significant scale dependence of the NLO re- 
sult. Similar discrepancies between the results of higher-order exact matrix element 
calculations and shower MC predictions have been reported in [Orr95, OrrQSa]. 

It should be pointed out that the transverse momentum of the top pair, directly 
related to the jet activity that accompanies a top production event, has important 
consequences on the determination of the experimental top detection efficiency, as 
well as a large impact on the reliability with which the mass of the top quark can be 
measured [Abe94, Abe95, Abachi95, Campagnari96] . The uncertainty related to the 
modelling of the gluon radiation emitted by initial and final state, in ti production 
and decay, is in fact one of the largest systematics present in the most recent top-mass 
determinations [Galtieri96, Grannis96]. 

Although the statistics of the data on top production are still low, the first distri- 
butions have been presented by CDF [Caner96]. The top-pair invariant-mass spec- 
trum is in good agreement with the theoretical expectations, and currently shows 
no evidence of an anomalous production source in the s-channel (fig. ^y). On the 
contrary, the pair transverse momentum is harder than the HERWIG MC predicts 



(fig. ^). As was pointed out above, this distribution is not, at the moment, predicted 
very accurately. It is therefore premature to conclude that the present discrepancy is 
physically significant, especially in view of the limited statistics. 

To summarize this section on top production, we conclude that nothing surprising 
has been discovered so far in the comparison of data with the predictions of NLO 
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Table 4: Average value of p!^ and fraction of the cross section with 
Py > 10 and 20 Ge V for different calculations. 
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Figure 40: CDF data on the invariant mass of top pairs, from W + 4 
jets (one of which tagged as a h), compared with the HERWIG MG 
expectations. 
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Figure 41: CDF data on the transverse momentum distribution of top 
pairs, from W + 4 jets (one of which tagged as a h), compared with the 
HERWIG MC expectations. 



perturbative QCD. As expected, the agreement between data and theory is aheady 
more sohd than in the case of bottom or charm production. Given the claimed accu- 
racy of the theoretical predictions, at least for the total cross section and for inclusive 
quantities, only additional statistics will however make these comparisons more com- 
pelling. Detailed features of the structure of the final state in top production, such 
as studies of the jet activity, have so far been probed only indirectly in the context 
of the mass measurements. It is likely that a lot will be learned in the future from 
these more detailed analyses, and that significant progress will be made in their MC 
modelling. 

5 Higher orders and resummation 



In this section we will deal with the problem of the resummation of logarithmi- 
cally enhanced effects in the vicinity of the threshold region in hard hadroproduction 
processes. Drell-Yan lepton-pair production has been in the past the best stud- 
ied example of this sort. The threshold region is reached when the invariant mass 
of the lepton pair approaches the total available energy. A large amount of theo- 
retical and phenomenological work has been done on this subject. The articles in 
refs. [Sterman87, Catani89] summarize the theoretical status of the subject, and 
also include references to the extensive literature in this field. For the Drell-Yan 
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process the resummation of soft gluons has been computed to NLO accuracy. Exten- 
sion of the NLO resummation formahsm to heavy-flavour production is under way 
[Kidonakis96, Contopanagos96]. 

Resummation formulae have also been used in estimating heavy-flavour produc- 
tion [Laenen92, Laenen94, Berger95, Berger96, Kidonakis95] . These works indicate 
the presence of very large higher-order corrections to heavy-flavour production at 
colliders, in particular in the gg production channel. Since top production at the 
Tevatron is dominated by the qq channel, the claimed size of the resummation ef- 
fects is about 15%, below the current experimental uncertainties on the cross section. 
Nevertheless, it is important to establish whether threshold resummation effects are 
really so large, a fact that would cast doubts on every QCD calculation of hadronic 
cross sections. Recently, in ref. [Catani96, Catani96a], it was argued that in fact re- 
summation effects are not as large as found in previous computations. In particular, 
it was pointed out there that large spurious terms are present in certain formulations 
of the resummation problem that are not justified by the threshold approximation. 
Furthermore, it was also shown that resummation may be formulated in such a way 
that these terms are absent. In the following we will present a brief review of this 
problem. 

5.1 What are soft-gluon effects 

Coloured particles emit soft gluons with high probability. Normally, the effect of 
soft-gluon emission is small (at least in inclusive quantities), since they only slightly 
affect the kinematics of a process. However, in a production process of high-mass ob- 
jects, when we approach the threshold, soft-gluon emission becomes important. Let 
us fix our attention on heavy-flavour production in hadronic collisions near threshold. 
The process is schematically depicted in fig. H2. The incoming protons make a big 
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Figure 42: Heavy-flavour production near ttireshold. 



effort in providing partons with a large fraction of the longitudinal momentum, thus 
going towards the very large x region of the structure functions. Under these circum- 
stances, even the small amount of energy wasted by soft-gluon radiation yields an 
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important suppression of the cross section. In the usual apphcation of the factoriza- 
tion theorem, part of the radiative corrections due to gluon radiation are included in 
the structure functions. Thus, depending upon the factorization scheme, the left-over 
suppression may be positive or negative. In the MS and DIS scheme the left-over 
is negative, so that the suppression effect appears instead as an enhancement of the 
cross section. 

In the case of heavy-flavour production, the perturbative expansion for the par- 
tonic cross section at 0(a|) (omitting obvious indices, such as the incoming parton 
types) has the structure 



a[s) = ao[s) 



l + Caslog'il-p) + 0{as\og{l-p),al) , (5.1 



where s is the square of the partonic centre-of-mass energy, p = Am? /s , cro(s) is the 
Born cross section 



a? 



^ol^J = ^/(°np) , (5.2) 

and as = as{rr?)- The reader can find in ref. [Nason88] explicit formulae for the 
functions f^^\p), as well as for the constant coefficient C. The precise value of C 
depends upon the type of incoming partons (i.e. quarks or gluons) and upon the 
factorization scheme. In both the MS and the DIS scheme it is positive, so that in 
the following we can think of it as being a positive constant. Resummation, according 
to refs. [Laenen92, Laenen94], gives 

a(^^^\s) = ^f^'\p)exp\Cas{s')\og\l-p)] , (5.3) 



m 

where s' is a scheme-dependent function of s that goes to zero as s-^4m^. We have 
s' = {1 — p)^m'^, where 1] = 1 in the DIS scheme, and 3/2 in the MS scheme. Formula 
( |5.3| ) is supposed to include all terms of order a^log"(l — p) with n > m in the 
exponent. Remember in fact that 

1 1 

"s(s') = — = — ^2 , , ,, 17 = tts (1 - asboV log(l - p) + . . .) , 

bo^ogj2 bo\og^ + bor]\og{l- p) 

(5.4) 

so that in the exponent there are terms with arbitrary powers of as, and a power of 

the logarithm that is always larger than the power of ag- The first subleading terms 

have the form a^ log'^(l — p). 

In order to get a physical cross section, the partonic cross section given above 
should be convoluted with parton luminosities: 

2 

a^'"'^ = ^ I C{r)f^'\p)exp[Cas{{l- prm')\og\l- p)]dT , (5.5) 

where, omitting obvious parton indices, 

C{r)=jFix,)FixM^,x,-r)dx,dx,. (5.6) 
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Here we can spot a problem in the resummation formula. When performing our 
integral over r, as p ^ 1 the argument of as in the exponential approaches zero. 
Before it actually hits the zero, it will hit a singularity in the running coupling as(s'), 
causing the integral to diverge. In order to avoid this divergence, a cutoff po was 
introduced in the literature [Laenen92, Laenen94, Appel88]. Observe that this cutoff 
has nothing to do with the standard factorization and renormalization scale /i. It is 
essentially a cutoff on soft-gluon radiation, imposed in order to avoid the blowing up 
of the running coupling associated with soft-gluon emission. 

The use of a cutoff seems an ad hoc procedure in this case. However, it can be 
justified to some extent. Suppose, for example, that we end up in a QCD calculation 
with a formula like 

G= f as{k^)G{k^)dk , (5.7) 

Jo 

where G{k'^) is a smooth function as k^ ^0. Integrals of this kind are often found, 
for example, in the computation of shape variables in jet physics. This expression 
is divergent as fc^— >-0, since at some point as approaches the Landau pole. The 
divergence can be handled by a cutoff po, which has to be large enough for as to be 
barely perturbative. For example, we may choose fiQ = 5A, a value around 2 GeV. 
We can then argue that 



G= [^ as{k^) GiP) dk + C^. (5 



In fact , the divergence of the 1-loop expression of as does not signal a real physical di- 
vergence. More likely, the point at which as becomes of order 1 signals the breakdown 
of perturbation theory. We therefore exclude this region, estimate its contribution by 
dimensional analysis, and obtain a power correction. A slightly more formal justifi- 
cation makes use of the concept of IR (infrared) renormalons. We expand eq. (|5.71 ) 
in powers of ag = as{Q'^), using 

1 cx) / u2\ i 

(5.9) 
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C^st[-^sbolog^'^ dk 

'^ POO 

= a^EK&o)^' / t^e-'/'dt, (5.10) 

where t = logQ'^/k'^. The integral can be performed, and one gets 



Goc2a,5^J!(2«A)^ (5.11) 

j=0 
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which is a divergent series, since a factorial grows faster than any power. This lack of 
convergence is in fact a general feature of the perturbative expansion in field theory. 
The perturbative expansion should be interpreted as an asymptotic one. The terms 
of the expansion ( ^.llj) decrease for moderate values of j. As j grows, the factorial 



takes over, the terms stop decreasing and begin to increase. This happens at the 
value of j at which the next term is equal to the current one 

ij + l)\{2asboy^'=j\i2asboy (5.12) 

or roughly 

^ (5.13) 



J mm — „ , • 

2a<j&o 

Asymptotic expansions are usually handled by summing their terms as long as they 
decrease. Of course, in this resummation prescription there is an ambiguity, which is 
of the order of the size of the first neglected term. In our case 

jmin! (2as&o)^"^'" ~ e^i°sJ-^l = e"^ ^ ^ , (5.14) 

which gives a power-suppressed ambiguity with the same power law that we found 
using the cutoff procedure. 

From the discussion given above, we would expect that the resummation formulae 
should include a cutoff of the order of a typical hadronic scale; varying the cutoff 
within a factor of order 1 should affect the cross section by terms of order A/Q. In 
fact, this is not the case. The cutoff has a dramatic effect on the cross section, as 
can be seen from figs. 2 and 3 of ref. [Laenen92]. For example, the uncertainty band 
obtained by varying the scale /iq between 0.2m and 0.3m, for top production in the 
gg channel, brings about a change in the cross section by a factor of 2, for a top 
mass between 100 and 200 GeV. These two values of ^o correspond to cutting off the 
gluon radiation at energies of the order of 20 to 30 GeV, therefore much larger than 
a typical hadronic scale. 

Other proposals for the resummation procedure have appeared in the literature. In 
refs. [Contopanagos94, Alvero95] a method was developed in the context of Drell-Yan 
pair production, which was applied to the heavy-flavour case in ref. [Berger95]. There, 
as can be seen from formula (116) and the subsequent discussion, unphysically large 
cutoffs are present, much larger than the typical hadronic scale one would expect. 

In the following section, we will show that the presence of large cutoffs and of 
large ambiguities in the resummation formula is not at all related to the blowing up 
of the coupling constant. In other words, there are other sources of factorial growth 
of the perturbative expansion for the resummation of soft gluons, and they largely 
dominate the factorial growth due to the running coupling. We will also show that 
these terms are spurious, and that soft-gluon resummation can be easily formulated 
in such a way that these terms are not present. 
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5.2 Problems with the aj-space resummation formula 

For definiteness, let us focus upon the resummation formula ( |5.5|) . We pointed 
out earlier that this formula is divergent when the argument of ag becomes too small 
and the coupling constant blows up. In fact, formula ( |5.5| ) is divergent, even for fixed 
coupling constant. At fixed coupling it can be written as 

„2 



^(rcs) ^ ^/rfr£(r)/W(p)expkClog2(l-/))" 

= —2!' dp 4^ (P/P) f^'\p) <^^P \o^sC\oi{l - p)] , (5.15) 

where p = Arri^/S, and the integral diverges as p—>-l, since the exponential 



exp(a log'(l - p)) = (1 - p)"i°gii-p) (5.16) 

grows faster than any inverse power of 1 — p as p — > 1 . This divergence can again 
be related to the factorial growth in the perturbative expansion. Expanding formula 
( |5.15| ) in powers of a^, we get 



E^(C«s)'7 dp-^^C{p/p)f^'\p)\og'\l-p). (5.17) 



^(res) ^ il£ y^ J_ 



It is now easy to see that, because of its singularity for p— ^1, the integral of log^^(l — p) 
grows like (2/c)! . Let us make here the simplifying assumption that /'•°^(p) ~ 9{l — p) . 
In a neighbourhood of the singularity, the integral behaves like 



/•l /"OO 

/ dp\og^\l-p)= dte-H^^ , {5.U 

Jl-e Jlogl/e 



where we have performed the substitution t = \ log(l — p)|. The above integral equals 
dt e-*P = / dt e-V' - / dte-V" = (2A;)! - / dte-'t^'' . (5.19) 

logl/e Jo Jo Jo 

Since 

r^ ^'^ dt e-H^'' < (log l/e)2'=+i (5.20) 

Jo 

we see that the contribution to the integral near the singularity is dominated by the 
term (2/c)! . The power expansion for the cross section is then 



^-' ^ E ^ (Cas)' ^ £ A:! (ACas)' . (5.21) 

fc=0 '^^ k=0 



00 
r(res) ^ 



The above formula is in fact appropriate for the case of the heavy-fiavour cross sec- 
tion with a lower cut on the invariant mass of the pair[5 As in the previous example, 

^^ In the case of the total cross section, using the known behaviour of /*^") near threshold, one 



would get an extra factor of 4/9 in front of C in eq. (5.21) 
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the factorial growth of formula ( 5.21|) will give rise to ambiguities in the resumma- 
tion of the perturbative expansion. These ambiguities are not, however, related to 
renormalons, since they occur also at fixed coupling constant. In fact, they are in 
general fractional powers of A/Q, where Q is the scale involved in the problem, and 
A is a typical hadronic scale. For example, as shown in ref. [Catani96a], in the case 
of heavy-flavour production through the gluon-fusion mechanism, at fixed invariant 
mass of the heavy-quark pair, the ambiguity would have the form (A/Q)°'^^, and it 
would thus be extremely relevant, even for very massive heavy-quark pairs. 

These large terms in the perturbative expansion are in fact spurious. They are an 
artefact of the x-space resummation procedure. This can be easily understood with 
the following argument. Exponentiation of the gluon emission is possible because, 
roughly speaking, each soft gluon is emitted independently. However, this indepen- 
dence is only approximate: the total momentum must be conserved. Momentum 
conservation, however, is a subleading effect in the soft resummation formula. Yet, 
its violation leads to factorially growing terms. These terms are subleading from the 
point of view of the logarithmic behaviour, but very important from the point of view 
of the factorial growth of the perturbative expansion. The presence of large factorial 
terms due to momentum non-conservation can be understood also by simple argu- 
ments. The emission of k gluons, where each gluon has a limit on its energy Ei < t], 
leads to a phase space that is larger by a factor of k\ than the case when the total 
energy of emission is bounded J^Ei < rj. Thus phase space alone provides a k\ factor. 
We can see in more detail the origin of the {2k)\ term by considering the formula 
for the partonic cross section with two emitted soft gluons, implementing momentum 
conservation. We have 



2 m^ 



log(l - z,] 

l-zi 



log(l - Z2) 

1-Z2 



f^>{p')6{p-ZiZ2p')dz,dZ2dp'. 



The leading-logarithmic term of the above integral is given by 



(5.22) 



n2«S 



a^'\p) = i-{2Casy^\og\l - p) f^'\p) + . . . 



(0)/ 



m^ 



(5.23) 



where terms with less than 4 powers of logarithms are neglected. We now see that 
the integral of the leading-logarithmic term of a^'^\p) in p has a large factor ^ 4! due 
to the integral of log^(l — p), while the integral of the full expression, eq. ( p. 22] ), gives 



dpa^^\p) = -{2Casf 



logfl 




dprKp) =0, (5.24) 



since the integrals of the soft emission factors vanish with the + prescription. In 
general, we see that if we take generic moments of a^^\p) (i.e. J p"^dpa^''\p)), the 
leading-log expression grows like (2^)! , while the full expression grows only geomet- 
rically with k. 
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The criticism described so far applies to the calculations of soft-gluon effects in 
heavy-flavour production given in refs. [Laenen92, Laenen94] and [Berger95]. As 
one may expect, the large factorial terms give rise to large corrections to the cross 
section. Since the terms of the perturbative expansion grow strongly with the order, 
they also give large uncertainties. In refs. [Laenen92, Laenen94], the presence of 
large uncertainties is in fact recognized. In ref. [Berger95] it is claimed that the 
uncertainties are small. Even there, however, an unphysically large cutoff is needed 
in order to make sense out of the resummation formulae. 

A second, more subtle problem with resummation prescriptions has to do with 
the presence of 1/Q corrections that arise from infrared renormalons. It was shown 
in ref. [Beneke95] that soft-gluon resummation does not yield the correct renormalon 
structure of the Drell-Yan cross section. More precisely, it was shown that in a 
fully calculable model the soft-gluon approximation yields a 1/Q correction, while 
in the exact result 1/Q effects are absent. This result suggests the absence of 1/Q 
corrections in Drell-Yan cross sections, an issue that is still much debated in the liter- 
ature [Korchemsky95 , Dokshitzer95, Akhoury95, Akhoury95a, Akhoury96, Beneke96, 
Korchemsky96]. We would like to remark, however, that even if 1/Q corrections were 
present in Drell-Yan and heavy-flavour production, they would only be of the order 
of 1% for top production at the Tevatron, and also that the result of ref. [Beneke95] 
shows that the soft-gluon approximation cannot describe them adequately. 

It is possible to formulate the resummation of soft gluons in such a way that the 
kinematic constraints are explicitly satisfied, and no factorial growth arises in the 
perturbative expansion. It is enough to formulate the resummation problem in the 
Mellin transform space^. 

In ref. [Catani96a] a resummation prescription using the Mellin space formula is 
given. It is shown there that there are, with this prescription, no factorially growing 
terms in the resummed perturbative expansion. Yet, soft effects are consistently in- 
cluded. This prescription is called "Minimal Prescription" (MP hereafter), because it 
does not introduce large terms that are not justified by the soft-gluon approximation. 

5.3 Phenomenological applications 

As stated earlier, when using the MP approach, one finds negligible resummation 
effects in most experimental configurations of interest. In figs. ^, ^ and ^, we plot 
the quantities 



"gg "qq Ogg + Oqq 



(5.25) 

\HH J 

.,^^ ..^^ ..^^ "^NLO 

Here S is equal to the MP-resummed hadronic cross section in which the terms of 
order q;| and a| have been subtracted, and unlo is the full hadronic NLO cross 
section. Thus, these plots show how important the impact of resummation is beyond 
the already computed NLO terms. The results for b production at the Tevatron 



^In fact, resummation formulae are usually derived in Mellin space. 
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Figure 43: Contribution of gluon resummation at order a|, and higher, 
relative to the NLO result, for the individual channels and for the total, 
for bottom production as a function of the CM energy in pp collisions. 
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Figure 44: Contribution of gluon resummation at order a| and higher, 
relative to the NLO result, for the individual subprocesses and for the 
total, as a function of the quark mass in pp collisions at 1.8 TeV. 
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at order a| and higher, relative to the NLO result, for the individual 
channels and for the total, as a function of the CM energy in pp colli- 
sions. 



can easily be inferred from fig. ^, since the qq component is negligible at Tevatron 
energies. 



For top production (figs. ^ and ^), we see that in most configurations of practical 
interest, the contribution of resummation is very small, being of the order of 1% at 
the Tevatron. A complete review of top-quark production at the Tevatron, based 
upon these findings, has been given in ref. [Catani96]. We also observe that, for top 
production at the LHC, soft-gluon resummation effects are equally negligible. Of 
course, in this last case, there are other corrections, not included here, that may 
need to be considered. Typically, since the values of x involved are small in this 
configuration, one may have to worry about the resummation of small-x logarithmic 
effects [Ellis90, Catani90, CataniQl, Collins91]. 

We see from the figures that in most experimental configurations of interest these 
effects are fully neghgible. One notable exception is b production at HERA-B, at 
yS = 39.2 GeV, where we find a 12% increase in the cross section. This correction is 
however well below the uncertainty due to higher-order radiative effects. For example, 
from the NLO calculation with the MRSA' [Martin94] parton densities and rrib = 
4.75 GeV, we get cXfej = 10.5^^'^ nb, a range obtained by varying the renormalization 
and factorization scales from mh/2 to 2mh. Thus the upper bound is 80% higher than 
the central value, to be compared with a 10% increase from the resummation effects. 
This result is much less dramatic than the results of ref. [Kidonakis95], obtained using 
the resummation prescription outlined in [Laenen92]. 
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6 Heavy-flavour production in e+e" collisions 

6.1 Preliminaries 

Heavy-flavour physics is a considerable part of tlie pliysics programme in e"'"e~ 
collisions, both at LEP and the SLC An entire chapter of this book is devoted 
to this topic [Kuhn97]. We shall limit ourselves here to some QCD aspects of the 
production dynamics. In particular we will discuss the most recent studies of the 
b fragmentation function, the quark-antiquark correlations, and the gluon-splitting 
production mechanism. 

We will not discuss specific topics related to the LEP2 programme. We would 
like, however, to point out that LEP2 offers the possibility of studying heavy-flavour 
production in photon-photon collisions. Relevant NLO computations have been per- 
formed for both the 77 [Drees93, Kraemer96] and the 67 cases [Laenen96]. Some 
phenomenological studies have already begun [Aurenche96]. A calculation of heavy- 
flavour production at high transverse momentum has also been presented [Cacciari96a]. 

6.2 Fragmentation function 

Theoretical predictions for fragmentation functions rely essentially upon two ingre- 
dients: perturbation theory, and a model for non-perturbative fragmentation effects. 
A comprehensive discussion of these topics is given in ref . [Nason92] . Leading-order 
formulae have been available for a long time [Azimov82, Azimov84], while the compu- 
tation of the NLO perturbative corrections is given in refs. [Mele90, Mele91]. Theoret- 
ical studies of the effects of non-perturbative physics are given in refs. [Colangelo92, 
Randall95, Nason96a]. Here we would like to call attention to the study of ref. 
[Colangelo92] , where a definite prediction was made for the 6- fragmentation func- 
tion at LEP for various values of Aqcd. In fig. ^ we report the predictions of 
ref. [Colangelo92] together with the data of ref. [Buskulic95] . Observe that the same 
figures were reported in ref. [Nason92], together with some unpublished L3 data, and 
at that time the agreement with theory seemed quite poor. Interestingly enough, the 
position of the peak is well reproduced by the QCD calculation. In general, a value of 
Aqcd between 200 and 300 MeV gives an adequate fit of the data. In ref. [Alexander96] 
a more detailed fit was performed, using the calculation of ref. [Colangelo92] . The 
results are shown in fig. |5^. There we can see that commonly used parametrizations 
fail to reproduce the slope of the fragmentation function, while the QCD result is 
quite consistent with the data. It also turns out that for the larger values of A the 

non-perturbative part of the initial condition is strongly peaked towards large values 

(5) 
of X. This means that, if we can assume that Aqcd > 0.25 GeV, the fragmentation 

function is largely perturbative, that is to say, is mostly made up of QCD evolu- 
tion and soft-gluon emission effects [Mele91]. The results of LEP measurements of 
the fragmentation functions are important also in view of possible applications to 
collider physics. The problems that are encountered in large transverse momentum 
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Figure 46: Theoretical 'prediction for ttie b fragmentation functions at 
LEP together with ALEPH data. Only the statistical error is shown. 



heavy-flavour production at hadron colliders may be partially due to a lack of proper 
treatment of the fragmentation effects. The study of the fragmentation function has 
concentrated in the past on the determination of its second moment (i.e. the average 
momentum of the B meson), as measured in e+e^ collisions. On the other hand, 
as pointed out in subsection |4.1| , h production at hadron colliders is also sensitive to 
higher moments of the fragmentation function. Since the new data from LEP exper- 
iments are now also sensitive to higher moments, it would be interesting to perform 
a new computation of the heavy-flavour production cross section at large transverse 
momentum, using fragmentation functions that flt the LEP data adequately. 



6.3 Heavy-quark production via gluon splitting 

The importance of the measurement of Ffe (the width of the Z^ into hh) has 
prompted several experimental and theoretical studies of 6-production characteris- 
tics (see, for instance, [Kuhn97]). One important issue is to understand how often hh 
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Figure 47: OPAL results for the h fragmentation function at LEP. 



or cc pairs are produced indirectly, via a gluon splitting mechanism. Several theoreti- 
cal studies are available on this topic [Nason92, Mueller86, Mangano92a, Seymour95]. 
Experimental studies on charm production via gluon splitting have been presented 
in refs. [Akers95, Akers95a, Hansper96], and a first measurement of g^bb has been 
given in [Branchini96]. The reported values are given in table §. In ref. [Seymour95] 
an explicit calculation of these quantities has been performed. Using these results^ 
we computed the charm and bottom multiplicities for different values of the masses 
and of Af^. We report them in table ^. As can be seen, the averaged experimental 
result of 2.38 ± 0.48% [Przysiezniak96] is consistent with the upper range of the theo- 
retical prediction, preferring lower values of the quark mass, and larger values of A^^. 



A similar problem, although much more severe, was mentioned in subsection [4.3.2 , 
in relation to charm jets in hadronic collisions. 

As reported in ref. [Seymour95], Monte Carlo models are in qualitative agreement 
with these results, although the spread of the values they obtain is somewhat larger 
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We thank M. Seymour for providing us with the relative FORTRAN code. 
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ng-,c-c (%) 


ng^bb (%) 


[Akers95a] 


2.27 ±0.28 ±0.41 




[Hansper96] 


2.65 ±0.74 ±0.51 




[Branchini96] 




0.22 ±0.10 ±0.08 


[Seymour95] 
Af = 150 MeV 
Af = 300 MeV 


1 qp;+0.48 
1 OCT +0.69 

l-85_0 44 


0.20 ±0.02 
0.26 ±0.03 



Table 5: Fraction of events containing g^cc and g^bb subprocesses 
in Z decays, as measured by the various collaborations, compared with 
theoretical predictions. The central values for the theoretical predictions 
are obtained with nic = 1-5 and nib = 4.75 GeV, the upper limits with 
rUc = 1.2 and nib = 4.5 GeV, and the lower limits with nic = 1.8 and 
nib = 5 GeV. 



than the theoretical error estimated by the direct calculation. In particular, one finds 
that while HERWIG [Marchesini92] and JETSET [Sjostrand87] agree quite well with 
the theoretical calculation, ARIADNE [Lonnblad92] is higher by roughly a factor of 
2, and thus is in better agreement with data. 



6.4 Correlations 



Correlations in the kinematical properties of the b and b also affect the systematics 
of the measurement of Tb- The most important effect is the correlation in the mo- 
mentum of the b and b quark, since detection efficiency is often strongly momentum- 
dependent. In ref. [Nason96] an explicit leading-logarithmic calculation of the mo- 
mentum correlation is given. The main result of ref. [Nason96] is the 0{as) formula 
for the double-inclusive distribution 



da 



dxidx2 



D{xi)D{x2) + —^ / dyidy2 9{yi ± 1/2 - 1) 7- 
37r Jo (1 



yl ± yl 



l/l)(l-y2) 




e{yi - xi) e{y2 - X2) - D{xi)D I—] - 9{y2 - X2) 

y2j yiy2 \y2 y2 



D{x2)-9{yi 

yi 



Xl, 



D{xi)D{x2) 



(6.1) 



where D{x) is the total heavy-flavour fragmentation function, and Xi and X2 are the 
energy fractions of the two heavy flavoured hadrons (x = 2E/yS). As an illustration. 
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Figure 48: Double inclusive cross section da/dxi dx2, plotted as a 
function of xi for several values of X2- 



we plot in fig. ^ the double-inclusive cross section da/dxi dx2 as a function of Xi for 
several values of X2- We use the Peterson parametrization given in eq. ( |2.9| ) with 
e = 0.04, which gives (x) = 0.70 and ag = 0.12. The positive momentum correlation 
is quite visible in fig. ^. As X2 increases, the peak of the distribution in xi also moves 
towards larger values. A particularly interesting quantity is the average momentum 
correlation, defined as 



(xiXa) - {x^ 



(6.2) 



This quantity is independent of the fragmentation function, and has therefore an 
expansion in the strong coupling with finite coefficients, starting at order as- In fact, 
one gets (in the limit rn^/S-^O) 



O.las + O{a 



(6.3) 



Since the correlation is so small, non-perturbative effects may be competitive with 
the perturbative ones. For example, if corrections of the order of Aqcd/v^ were 
present, this would be the case at LEP energies. In ref. [Nason96a] it was shown that 
in the renormalon approach power corrections to correlations are suppressed at the 
level of Aqqj-,/S' at least. This result, being based upon the renormalon approach, is 
not fully conclusive. It is however encouraging, since it suggests that the perturbative 
calculation should be very reliable at LEP energies. 



84 Heavy- Quark Production 



7 Conclusions and outlook 

We collect here the main conclusions of the various sections, and present our 
outlook on future progress in the field. 

The experimental results on total cross sections for charm production at fixed 
target are in good agreement with NLO QCD. In the case of hadroproduction the ex- 
perimental accuracy is far superior to the theoretical accuracy, and might ultimately 
help in better pinning down the values of the theoretical parameters. For example, 
current data in pion- and proton-induced reactions strongly disfavour a charm mass 
value as large as 1.8 GeV. Similar conclusions could not be drawn from a study of 
the charm photoproduction data. This is because, although the general precision 
of the data is better than the overall theoretical uncertainty, data from different 
experiments have a large spread, not obviously compatible with the statistical and 
systematic uncertainties quoted by the single experiments. Within the large theo- 
retical and experimental uncertainties, the agreement with theoretical expectations 
is nevertheless good, and consistent with a charm mass value of 1.5 GeV, and with 
parameters such as A^^ and PDF sets favoured by the hadroproduction studies. 

Bottom production in fixed-target experiments is affected by the low statistics 
available. Even within the statistical uncertainties, there are however conflicting 
results in the present data. The agreement with theory, afflicted by equally large 
uncertainties, is therefore at this time not particularly enlightening. These conflicts 
should be resolved by the current and next generation of high- sensitivity experiments 
(at FNAL and HERA-B). 

While the theoretical description of the total production rates is expected to be 
insensitive to possible non-perturbative aspects of the production mechanisms, differ- 
ential distributions for charmed hadrons produced in flxed-target experiments provide 
an important probe of these phenomena. The overall conclusions of the studies we 
presented can be summarized as follows. The inclusive pj. spectra predicted by NLO 
QCD for charm photoproduction are too hard to agree with the data; they show 
clear evidence of a non-perturbative fragmentation mechanism that slows down the 
charmed hadrons. If this fragmentation is then applied to charm hadroproduction, 
the spectra turn out to be far too soft. The apparent inconsistency can be solved 
by invoking the effect of an intrinsic non-perturbative /cj. kick of the partons inside 
the hadron. This kr kick makes the pr spectrum harder, an effect that is stronger in 
hadroproduction than in photoproduction because of the presence of two hadrons in- 
stead of one in the initial state. The need for a substantial /c^ kick, of the order of 0.7 
to 1 GeV, is corroborated by the study of azimuthal correlations of charm pairs, both 
in photo- and in hadroproduction. There remain indications that a non-perturbative 
fragmentation slightly softer than the standard Peterson parametrization would pro- 
vide a better agreement with data. Nevertheless the overall picture is that once these 
two main effects are accounted for, the qualitative and quantitative features of most 
of the available distributions are properly described by the theory, independently of 
the nature of the beam. 
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Additional non-perturbative effects, such as colour-drag from the beam, should be 
invoked to properly describe the Xp distribution at large Xp and D/D asymmetries. 
The modelling of colour-drag effects contained in the available Monte Carlo programs 
can be properly tuned to provide an acceptable description of the data. 

Heavy-quark photo- and electroproduction studies at HERA are still in an early 
stage, with much more statistics waiting to become available in the near future. The 
current data on charm total cross sections agree with the NLO QCD calculations for 
choices of parameters consistent with low-energy fixed-target results. However, the 
current HERA experimental accuracy and the spread of the low-energy photoproduc- 
tion data are not sufficient as yet to improve our knowledge of the gluon densities 
of either the proton or the photon. Differential distributions in p^ and in rapidity 
have also been presented, and a general agreement with the theory was found. Some 
discrepancies, as found for example in the rapidity distributions, are not statistically 
compelling as yet. 

Progress in the study of the gluon density of the proton will come when higher 
statistics, and in particular large samples of charm pairs, will be available. Improved 
efficiencies at large rapidity and higher statistics will likewise help the understanding 
of the gluon structure of the photon. 

The study of bottom production at the Tevatron shows that there is good agree- 
ment between the shape of the 6-quark p^ distribution predicted by NLO QCD and 
that observed in the data for central rapidities. A similar conclusion can be reached 
in the case of the shape of the azimuthal correlations between b quarks. Although 
the data rate is higher by a factor of approximately 2 relative to the default choice 
of theoretical parameters, acceptable choices of A^'^ and of renormalization and fac- 
torization scales bring the theory in perfect agreement with the data of UAl and 
DO, and within 30% of the CDF measurements. It is encouraging that studies of 
higher-order logarithmic corrections favour the choice of low values for /iR and ^p. 
The CDF measurements at 630 and 1800 GeV indicate that theory correctly predicts 
the scaling of the differential 6-quark p^ distribution between 630 and 1800 GeV, a 
fact that had often been questioned in the past and now finds strong support. 

More theoretical studies should be devoted to the understanding of the non- 
perturbative fragmentation function for heavy quarks. We showed that the inter- 
pretation of the experimental data on the p^ spectrum depends very strongly on the 
assumed shape of the fragmentation function. Measurements of the 6-quark content of 
high-i?y jets, a quantity that is rather independent of the fragmentation properties of 
the b quark, indicate a good agreement between data and theory. This suggests that 
residual discrepancies in the comparison of theory and data for the p^ spectra, and 
possibly for the azimuthal correlations, could be resolved with a better understanding 
of the fragmentation phenomena. 

No surprise has so far come up in the comparison of the data on top production 
with the predictions of NLO perturbative QCD. Given the claimed accuracy of the 
theoretical predictions, which in the case of the total cross section is of the order 
of 10%, only additional statistics will make these comparisons more compelling. We 
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remark here that, on the theoretical side, soft-gluon resummation effects have been 
found to be much less important than previously thought. Furthermore, the ag de- 
pendence of the top cross section is quite small, because of a compensation mechanism 
due to structure-function evolution effects. Thus, the theoretical error on the cross 
section is quite credible in this case. Detailed features of the structure of the final 
state in top production, such as studies of the jet activity, have so far been probed 
only indirectly in the context of the mass measurements. It is likely that a lot will be 
learned in the future from these more detailed analyses, and that significant progress 
will be made in their MC modelling. 

Considerable progress has been made in the study of heavy-fiavour production in 
e^e^ annihilation. The measurement of the B fragmentation function at LEP has 
reached a remarkable level of accuracy. Some details of the production mechanism, 
particularly important for the determination of F^ have received special attention. 
Thus, the heavy-fiavour production by the gluon-splitting mechanism has been stud- 
ied both experimentally and theoretically, and analytical calculations of the corre- 
lations of the heavy-fiavour pair have been performed. Besides their intrinsic value, 
e^e~ studies could clarify several aspects of the production mechanism, which are 
important also at hadron colliders. 
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